Source: THE PROPERTIES OF GASES AND LIQUIDS 


CHAPTER 9 

VISCOSITY 


9-1 Scope 


The first part of this chapter deals with the viscosity of gases and the second with 
the viscosity of liquids. In each part, methods are recommended for: (1) correlating 
viscosities with temperature; (2) estimating viscosities when no experimental data 
are available; (3) estimating the effect of pressure on viscosity; and (4) estimating 
the viscosities of mixtures. The molecular theory of viscosity is considered briefly. 


9-2 Definitions of Units of Viscosity 


If a shearing stress is applied to any portion of a confined fluid, the fluid will move 
with a velocity gradient with its maximum velocity at the point where the stress is 
applied. If the local shear stress per unit area at any point is divided by the velocity 
gradient, the ratio obtained is defined as the viscosity of the medium. Thus, vis¬ 
cosity is a measure of the internal fluid friction, which tends to oppose any dynamic 
change in the fluid motion. An applied shearing force will result in a large velocity 
gradient at low viscosity. Increased viscosity causes each fluid layer to exert a larger 
frictional drag on adjacent layers which in turn decreases the velocity gradient. 

It is to be noted that viscosity differs in one important respect from the properties 
discussed previously in this book; namely, viscosity can only be measured in a 
nonequilibrium experiment. This is unlike density which can be found in a static 
apparatus and so is an equilibrium property. On the microscale, however, both 
properties reflect the effects of molecular motion and interaction. Thus, even though 
viscosity is ordinarily referred to as a nonequilibrium property; it is, like density, 
a function of the thermodynamic state of the fluid; in fact, it may even be used to 
define the state of the material. Brule and Starling (1984) have emphasized the 
desirability of using both viscosity and thermodynamic data to characterize complex 
fluids and to develop correlations. This discussion is limited to Newtonian fluids, 
i.e., fluids in which the viscosity, as defined, is independent of either the magnitude 
of the shearing stress or velocity gradient (rate of shear). For polymer solutions 
which are non-Newtonian, the reader is referred to Ferry (1980) or Larson (1999). 

The mechanisms and molecular theory of gas viscosity have been reasonably 
well clarified by nonequlibrium statistical mechanics and the kinetic theory of gases 
(Millat, et al., 1996), but the theory of liquid viscosity is less well developed. Brief 
summaries of both theories will be presented. 


9.1 
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Since viscosity is defined as a shearing stress per unit area divided by a velocity 
gradient, it should have the dimensions of (force) (time)/(length) 2 or mass/(length) 
(time). Both dimensional groups are used, although for most scientific work, vis¬ 
cosities are expressed in poises, centipoises, micropoises, etc. A poise (P) denotes 
a viscosity of 0.1 N ■ s/m 2 and 1.0 cP = 0.01 P. The following conversion factors 
apply to viscosity units: 

1 P = 100 cP = 1.000 X 10 6 p,P = 0.1 N s/m 2 = 1 g/(cm ■ s) = 0.1 Pa • s 
= 6.72 X 10- 2 lb-mass/(ft • s)= 242 lb-mass/(ft • h) 

1 cP = 1 mPa ■ s 

The kinematic viscosity is the ratio of the viscosity to the density. With viscosity 
in poises and the density of grams per cubic centimeter, the unit of kinematic 
viscosity is the stoke, with the units square centimeters per second. In the SI system 
of units, viscosities are expressed in N ■ s/m 2 (or Pa ■ s) and kinematic viscosities 
in either m 2 /s or cm 2 /s. 


9-3 Theory of Gas Transport Properties 


The theory of gas transport properties is simply stated, but it is quite complex to 
express in equations that can be used directly to calculate viscosities. In simple 
terms, when a gas undergoes a shearing stress so that there is some bulk motion, 
the molecules at any one point have the bulk velocity vector added to their own 
random velocity vector. Molecular collisions cause an interchange of momentum 
throughout the fluid, and this bulk motion velocity (or momentum) becomes dis¬ 
tributed. Near the source of the applied stress, the bulk velocity vector is high, but 
as the molecules move away from the source, they are “slowed down” due to 
random molecular collisions. This random, molecular momentum interchange is the 
predominant cause of gaseous viscosity. 


Elementary Kinetic Theory 

If the gas is modeled in the simplest manner, it is possible to show the general 
relations among viscosity, temperature, pressure, and molecular size. More rigorous 
treatments yield similar relations but with important correction factors. The ele¬ 
mentary gas model assumes all molecules to be nonattracting rigid spheres of di¬ 
ameter cr (with mass m) moving randomly at a mean velocity v. The molar density 
is n molecules in a unit volume while the mass density is the mass in a unit volume. 
Molecules move in the gas and collide transferring momentum in a velocity gradient 
and energy in a temperature gradient. The motion also transfers molecular species 
in a concentration gradient. The net flux of momentum, energy, or component mass 
between two layers is assumed proportional to the momentum, energy, or mass 
density gradient, i.e., 


Flux °c — 7 - (9-3.1) 

dz 

where the density p' decreases in the +z direction and p' may be p, (mass density), 
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nmv , (momentum density), or CjiT (energy density). The coefficient of propor¬ 
tionality for all these fluxes is given by elementary kinetic theory as vLI 3, where 
v is the average molecular speed and L is the mean free path. 

Equation (9-3.1) is also used to define the transport coefficients of diffusivity D, 
viscosity 17 , and thermal conductivity A; that is. 


Mass flux = 


„ dn, vL dp, 

-Dm — = — — 

dz 3 dz 


dv. 


vL 


dv„ 


Momentum flux = — n —r"- =-— mn -p- 

3 dz 


dz 


Energy flux = 


dT vL „ dT 

“A — = - — C u n — 
dz 3 dz 


(9-3.2) 

(9-3.3) 

(9-3.4) 


Equations (9-3.2) to (9-3.4) define the transport coefficients I), 17 , and A. If the 
average speed is proportional to ( RT/M) 1 ' 2 and the mean free path to ( ncr 2 ) -1 , 


„ vL , , r 3/2 

D = - = (const) —— 

mpvL T 1/2 M m 

V = —T~ = (const)- i — 

3 (T 

vLc v n T 1 ' 2 

A = _ = ( C ons t) __ 


(9-3.5) 

(9-3.6) 

(9-3.7) 


The constant multipliers in Eqs. (9-3.5) to (9.3-7) are different in each case; the 
interesting fact to note from these results is the dependency of the various transfer 
coefficients on T, P, M, and <r. A similar treatment for rigid, nonattracting spheres 
having a Maxwellian velocity distribution yields the same final equations but with 
slightly different numerical constants. 

The viscosity relation [Eq. (9-3.6)] for a rigid, non-attracting sphere model is 
(see page 14 of Hirschfelder, et al., 1954) 

(MT) 112 

7 ) = 26.69 2 -(9-3.8) 

<x 


where 77 = viscosity, pP 

M = molecular weight, g/mol 
T = temperature, K 
cr = hard-sphere diameter, A 

Analogous equations for A and D are given in Chaps. 10 and 11. 


Effect of Intermolecular Forces 

If the molecules attract or repel one another by virtue of intermolecular forces, the 
theory of Chapman and Enskog is normally employed (Chapman and Cowling, 
1939; Hirschfelder, et al., 1954). There are four important assumptions in this de¬ 
velopment: ( 1 ) the gas is sufficiently dilute for only binary collisions to occur (ideal 
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gas); (2) the motion of the molecules during a collision can be described by classical 
mechanics; (3) only elastic collisions occur, and (4) the intermolecular forces act 
only between fixed centers of the molecules; i.e., the intermolecular potential func¬ 
tion is spherically symmetric. With these restrictions, it would appear that the re¬ 
sulting theory should be applicable only to low-pressure, high-temperature mona¬ 
tomic gases. The pressure and temperature restrictions are valid for polyatomic 
gases and, except for thermal conductivity (See Chap. 10) are adequate for most 
modeling purposes. 

The Chapman-Enskog treatment develops integral relations for the transport 
properties when the interactions between colliding molecules are described by a 
potential energy function <//(r). The equations require complex numerical solution 
for each choice of intermolecular potential model. In general terms, the first-order 
solution for viscosity can be written 


„ = (26 - 6 ^ r)1/2 (9-3.9) 

where the temperature dependence of the collision integral, IT,, is different for each 
<//(r) and all symbols and units are as defined in Eq. (9-3.8). fl v is unity if the 
molecules do not attract each other. Corrections can be found in Chapman and 
Cowling (1939) and Hirschfelder, et al. (1954). The use of fl„ from the Lennard- 
Jones (12-6) potential function is illustrated in Sec. 9-4. 


9-4 ESTIMATION OF LOW-PRESSURE GAS 
VISCOSITY 


Essentially all gas viscosity estimation techniques are based on either the Chapman- 
Enskog theory or the law of corresponding states. Both approaches are discussed 
below, and recommendations are presented at the end of the section. Experimental 
values of low-pressure gas viscosities are compiled in Landolt-Bornstein (1955), 
Stephan and Lucas (1979), and Vargaftik, et al. (1996). Literature references for a 
number of substances along with equations with which to calculate gas viscosities 
based on critically evaluated data may be found in Daubert, et al. (1997). Gas phase 
viscosity information can also be found in Dean (1999), Lide (1999), Perry and 
Green (1997), and Yaws (1995, 1995a). This information should be used with cau¬ 
tion in those cases where constants in equations have been determined from esti¬ 
mated rather than experimental viscosities. 


Theoretical Approach 

The first-order Chapman-Enskog viscosity equation was given as Eq. (9-3.9). To 
use this relation to estimate viscosities, the collision diameter cr and the collision 
integral must be found. In the derivation of Eq. (9-3.9), 11,, is obtained as a 
function of a dimensionless temperature T* which depends upon the intermolecular 
potential chosen. For any potential curve, the dimensionless temperature T* is re¬ 
lated to £ by 


kT 

T* — _ 


£ 


(9-4.1) 
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where k is Boltzmann’s constant and e is the minimum of the pair-potential energy. 
The working equation for 17 must have as many parameters as were used to define 
the original i/r(r) relation. While many potential models have been proposed 
(Hirschfelder, et ah, 1954), the Lennard-Jones 12-6 was the first and has most often 
been applied for ideal gas viscosity. 


i/r(r) = 4e 




(9-4.2) 


In Eq. (9-4.2), cr is like a molecular diameter and is the value of r that causes 
i l>(r) to be zero. With this potential, the collision integral has been determined by 
a number of investigators (Barker, et al., 1964; Hirschfelder, et al., 1954; Itean, et 
al., 1961; Klein and Smith, 1968; Monchick and Mason, 1961; and O’Connell and 
Prausnitz, 1965). Neufeld, et al. (1972) proposed an empirical equation which is 
convenient for computer application: 

Xl„ = [A(T*y B ] + C[exp(-DT *)J + Etexpt-FT*)] (9-4.3) 

where T* = kT/e, A = 1.16145, B = 0.14874, C = 0.52487, D = 0.77320, E = 
2.16178, and F = 2.43787. Equation (9-4.3) is applicable from 0.3 < T* < 100 
with an average deviation of only 0.064%. A graph of log 1 \ as a function of log 
T* is shown in Fig. 9-1. 



- 0.5 - 0.4 - 0.3 - 0.2 - 0.1 0.0 0.1 


log r* 

FIGURE 9-1 Effect of temperature on the Lennard-Jones viscosity 
collision integral. 
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With values of Cl v as a function of T*, a number of investigators have used Eq. 
(9-3.9) and regressed experimental viscosity-temperature data to find the best values 
of elk and cr for many substances. Appendix B lists a number of such sets as 
reported by Svehla (1962). It should be noted, however, that there appears also to 
be a number of other quite satisfactory sets of elk and cr for any given compound. 
For example, with n-butane, Svehla suggested e Ik = 513.4 K, cr = 4.730 A, 
whereas Flynn and Thodos (1962) recommend s/k = 208 K and cr = 5.869 A. 

Both sets, when used to calculate viscosities, yield almost exactly the same 
values of viscosity as shown in Fig. 9-2. This interesting paradox has been resolved 
by Reichenberg (1971), who suggested that log f l D is essentially a linear function 
of log T* (see Fig. 9-1). 


(1„ = a ( T*f (9-4.4) 

Kim and Ross (1967) do, in fact, propose that: 

n u = 1.604 (T*)-°- 5 (9-4.5) 

where 0.4 < T* < 1.4. They note a maximum error of only 0.7%. Substitution of 
Eq. (9-4.5) into Eq. (9-3.9) leads to 



270 300 350 400 


Temperature, K 

FIGURE 9-2 Comparison of calculated and exper¬ 
imental low-pressure gas viscosity of n-butane with 
Eq. (9-3.9) and the Lennard-Jones Potential: — Flynn 
and Thodos (1962) with cr = 5.869 A and elk = 208 

K. _Svehla (1962) with a = 4.730 A and elk 

= 513.4 K. o ( Titani, 1929 ); V (Wobster and Mueller, 
1941). 
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16.64(M) 1/2 r 
(. s/k) 1 , 2 o - 2 


(9-4.6) 


where the units are the same as in Eq. (9-3.9). Here the parameters a and slk are 
combined as a single term (E/k) 1,2 cr 2 . There is then no way of delineating individual 
values of e Ik and a by using experimental viscosity data, at least over the range 
where Eq. (9-4.5) applies. Equation (9-4.6) suggests that when 7) < 1.4, low- 
pressure gas viscosities are essentially proportional to the absolute temperature. 

The conclusion to be drawn from this discussion is that Eq. (9-3.9) can be used 
to calculate gas viscosity, although the chosen set of e Ik and a may have little 
relation to molecular properties. There will be an infinite number of acceptable sets 
as long as the temperature range is not too broad, e.g., if one limits the estimation 
to the range of reduced temperatures from about 0.3 to 1.2. In using published 
values of s/k and cr for a fluid of interest, the two values from the same set must 
be used—never e Ik from one set and cr from another. 

The difficulty in obtaining a priori meaningful values of s/k and cr has led most 
authors to specify rules which relate e Ik and a to macroscopic parameters such as 
the critical constants. One such method is shown below. 


Method of Chung, et al. (1984, 1988) 


These authors have employed Eq. (9-3.9) with 

e T c 
k 1.2593 


(9-4.7) 


o- = 0.809V,! 73 


(9-4.8) 


where e Ik and T c are in kelvins, cr is in angstroms, and V c is in cm 3 /mol. Then, 
using Eqs. (9-4.1) and (9-4.7), 


T* = 1.25937). 


(9-4.9) 


fl„ in Eq. (9-3.9) is found from Eq. (9-4.3) with T* defined by Eq. (9-4.9). Chung, 
et al. also multiply the right-hand side of Eq. (9-3.9) by a factor F c to account for 
molecular shapes and polarities of dilute gases. Their final result may be expressed 
as: 


V 


40.785 


F c (MT) m 

v™n D 


(9-4.10) 


where -p 
M 
T 

Vc 

F. 


viscosity, /j,P 
molecular weight, g/mol 
temperature, K 
critical volume, cm 3 /mol 

viscosity collision integral from Eq. (9-4.3) and T* = 1.25937) 

1 - 0.2756&J + 0.059035 + k (9-4.11) 


In Eq. (9-4.11), co is the acentric factor (See Chap. 2) and k is a special correction 
for highly polar substances such as alcohols and acids. Values of k for a few such 
materials are shown in Table 9-1. Chung, et al. (1984) suggest that for other al- 
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TABLE 9-1 The Association Factor k in Eq. 
(9-4.11) (Chung, et al., 1988) 


Compound 

K 

Compound 

K 

Methanol 

0.215 

n-Pentanol 

0.122 

Ethanol 

0.175 

n-Hexanol 

0.114 

n-Propanol 

0.143 

n-Heptanol 

0.109 

/-Propanol 

0.143 

Acetic Acid 

0.0916 

n-Butanol 

/-Butanol 

0.132 

0.132 

Water 

0.076 


cohols not shown in Table 9-1, k = 0.0682 + 4.704[(number of -OH groups)/ 
molecular weight]. The term /ju r is a dimensionless dipole moment. (See discussion 
in Chap. 2 and also under Eq. (9-4.17) for techniques to nondimensionalize a dipole 
moment.) When V c is in cm 3 /mole, T c is in kelvins, and /jl is in debyes, 

^ 13L3 (W (9 ' 4 - 12) 


Example 9-1 Estimate the viscosity of sulfur dioxide gas at atmospheric pressure 
and 300°C by using the Chung, et al. method. The experimental viscosity is 246 /rP 
(Landolt-Bomstein, 1955). 


solution From Appendix A, T c = 430.8 K, V c = 122 cm 3 /mole, M = 64.065, and 
the dipole moment is 1.6 debyes. From Eq. (2-3.3), to = 0.257. Assume k is negligible. 
From Eq. (9-4.12), 


(131.3)(1.6) 

[(122)(430.8)] 1/2 


0.916 


and with Eq. (9-4.11), 

F c = 1 — (0.27561(0.257) + (0.059035)(0.916) 4 = 0.971 
300 + 273 

r * = L2593 ^r = 1 - 675 


Then, with Eq. (9-4.3), fl v = 1.256. The viscosity is determined from Eq. (9-4.10). 

1(64.065X300 + 273)] 1 ' 2 

r, = (40.785X0.971) (122 ( )2/3( ,. 256) = 245.5 m P 

Error = 245 ' 5 ~ 246 x 100 = -0.2% 

246 


Experimental viscosities and those estimated by the Chung, et al. method are 
shown in Table 9-2. The critical properties used to prepare Table 9-2 differed 
slightly from those in Appendix A. The average absolute error was about 1.9%. 
This agrees well with the more extensive comparison by Chung (1980), who found 
an average absolute error of about 1.5%. 
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Corresponding States Methods 


From an equation such as (9-3.9), if one associates cr 3 with V c [as in Eq. (9-4.8)] 
and assumes V c is proportional to RT C /P C , a dimensionless viscosity can be defined: 


Vr = €v = f(T r ) 

= ~ (RT c )(N a ) 2 ] 1 

M 2 P 4 C 


(9-4.13) 

(9-4.14) 


In SI units, if R = 8314 J/(kmol • K) and N a (Avogadro’s number) = 6.023 X 10 26 
(kmolr 1 and with T c in kelvins, M in kg/kmol, and P c in N/m 2 , £ has the units of 
m 2 /(N • s) or inverse viscosity. In more convenient units, 

«- 0176 (thij ” (9 - 415) 


where f — reduced, inverse viscosity, (/xP) ', T c is kelvins, M is in g/mol, and P c 
is in bars. 

Equation (9-4.13) has been recommended by several authors (Flynn and Thodos, 
1961; Golubev, 1959; Malek and Stiel, 1972; Mathur and Thodos, 1963; Trautz, 
1931; and Yoon and Thodos, 1970). The specific form suggested by Lucas (Lucas, 
1980; Lucas, 1983; and Lucas, 1984a) is illustrated below. 


tj£ = [0.807 T° r 6ls - 0.357 exp(-0.449r,.) 

+ 0.340 exp(—4.058 T r ) + 0.018] F° P F° Q (9-4.16) 

£ is defined by Eq. (9-4.15), 17 is in /jlP, /', is the reduced temperature, and F° P and 
F° q are correction factors to account for polarity or quantum effects. To obtain F P , 
a reduced dipole moment is required. Lucas defines this quantity somewhat differ¬ 
ently than did Chung, et al. in Eq. (9-4.12), i.e., 

u?P 

IL r = 52.46 ^ (9-4.17) 


where /jl is in debyes, P c is in bars, and T c is in kelvins. Then F° P values are found 
as: 

F° P = 1 0 < n r < 0.022 

F° P = 1 + 30.55(0.292 - Z c ) 112 0.022 </«.,< 0.075 

F° = 1 + 30.55(0.292 - Z c ) 112 |0.96 + 0.1(7; - 0.7)| 0.075 < ix r 

(9-4.18) 

The factor F° Q is used only for the quantum gases He, H 2 , and D 2 . 

F° q = 1.22(2°' 15 {1 + 0.00385[(T r - 12) 2 ] 1/M sign (T, - 12)} (9-4.19) 

where Q — 1.38 (He), Q = 0.76 (H 2 ), Q = 0.52 (D 2 ). Sign( ) indicates that one 

should use +1 or — 1 depending on whether the value of the argument ( ) is greater 
than or less than zero. 
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TABLE 9-2 Comparison between Calculated and Experimental Low-Pressure Gas 
Viscosities 


Percent errort 


Chung 
et al., 


Compound 

T, °C 

Experimental 
value, /j.Pt 

Eq. 

(9-4.10) 

Lucas, 

Eq. (9-4.16) 

Reichenberg, 
Eq. (9-4.21) 

Acetic acid 

150 

118 

3.4 

— 

1.7 


250 

151 

0.2 

— 

2.4 

Acetylene 

30 

102 

0.6 

0.6 

2.4 


101 

126 

-0.6 

-0.8 

0.7 


200 

155 

-0.5 

-0.6 

0.7 

Ammonia 

37 

106 

2.2 

2.3 

— 


147 

146 

0.4 

0.2 

— 


267 

189 

-2.1 

-1.9 

— 

Benzene 

28 

73.2 

1.0 

3.2 

4.4 


100 

92.5 

0.2 

1.3 

2.1 


200 

117 

0.5 

1.6 

1.8 

Bromotrifluoromethane 

17 

145 

7.6 

8.3 

-0.3 


97 

183 

8.2 

8.6 

-0.5 

Isobutane 

37 

79 

-1.6 

1.6 

1.6 


155 

105 

2.0 

4.7 

4.0 


287 

132 

3.9 

6.7 

5.6 

n-butane 

7 

74 

-5.8 

-4.1 

-5.6 


127 

101 

-1.0 

0.1 

-2.3 


267 

132 

0.6 

1.5 

-1.4 

1-Butene 

20 

76.1 

-1.0 

1.1 

1.8 


60 

86.3 

-0.5 

1.4 

1.7 


120 

102 

-0.6 

0.9 

0.9 

Carbon dioxide 

37 

154 

-0.6 

1.6 

— 


127 

194 

-0.1 

2.2 

— 


327 

272 

0.5 

3.0 

— 

Carbon disulfide 

30 

100 

0 

5.9 

— 


98.2 

125 

-1.3 

3.5 

— 


200 

161 

-1.9 

2.2 

— 

Carbon tetrachloride 

125 

133 

0.8 

1.9 

-2.8 


200 

156 

2.4 

3.1 

-2.1 


300 

186 

3.8 

4.1 

-1.4 

Chlorine 

20 

133 

2.2 

3.5 

— 


100 

168 

3.3 

4.1 

— 


200 

209 

4.5 

5.1 

— 

Chloroform 

20 

100 

-3.8 

3.4 

-1.0 


100 

125 

-1.3 

5.4 

0.6 


300 

191 

-0.9 

5.1 

-0.6 

Cyclohexane 

35 

72.3 

-2.4 

-0.1 

-4.2 


100 

87.3 

-1.4 

0.3 

-4.1 


300 

129 

2.6 

3.6 

-1.9 

Dimethyl ether 

20 

90.9 

-5.2 

-1.0 

2.0 


100 

117 

-6.1 

-2.4 

0.2 
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TABLE 9-2 Comparison between Calculated and Experimental Low-Pressure Gas 
Viscosities ( Continued ) 


Compound 

/', °C 

Experimental 
value, /jlP$ 

Chung 
et al., 

Eq. 

(9-4.10) 

Percent errort 

Lucas, Eq. 
(9-4.16) 

Reichenberg, 
Eq. (9-4.21) 

Ethane 

47 

100 

0.1 

0.7 

-3.2 


117 

120 

0.2 

0.8 

-3.3 


247 

156 

-1.0 

-0.3 

-4.7 

Ethyl acetate 

125 

101 

-2.6 

9.0 

-1.5 


200 

120 

-2.4 

8.8 

-2.0 


300 

146 

-3.3 

7.4 

-3.5 

Ethanol 

110 

111 

-0.5 

-2.5 

0.6 


197 

137 

-0.8 

1.6 

-0.4 


267 

156 

-0.3 

2.0 

-0.1 

Diethyl ether 

125 

99.1 

-0.4 

0.2 

0.7 


200 

118 

-0.8 

-0.5 

-0.3 


300 

141 

-0.6 

-0.4 

-0.5 

n-Hexane 

107 

81 

-1.1 

0.9 

0.3 


267 

116 

-2.1 

-0.7 

-2.1 

Methane 

-13 

98 

-0.7 

-0.5 

— 


147 

147 

0 

-0.9 

— 

Methyl acetate 

125 

108 

0 

10 

1.7 


200 

157 

-1.6 

7.6 

-1.0 

Methanol 

67 

112 

-0.4 

-0.9 

1.1 


127 

132 

-0.3 

-1.1 

1.2 


277 

181 

-0.3 

-1.7 

0.9 

Methyl chloride 

50 

119 

3.9 

0.8 

-0.6 


130 

147 

5.0 

1.4 

1.9 

Nitrogen 

27 

178 

0.3 

0.3 

— 


227 

258 

0.3 

0.4 

— 

Isopropanol 

157 

113 

0.3 

5.9 

3.5 


257 

139 

0.1 

5.4 

2.7 

Propylene 

17 

83 

2.4 

2.4 

3.0 


127 

115 

1.5 

0.9 

0.8 


307 

160 

1.7 

1.2 

0.7 

Sulfur dioxide 

10 

120 

2.8 

5.5 

— 


100 

163 

0.3 

2.2 

— 


300 

246 

-0.2 

1.5 

— 


700 

376 

1.5 

3.4 

— 

Toluene 

60 

78.9 

-5.2 

-3.7 

-2.3 


250 

123 

-3.5 

-3.1 

-2.5 

Average absolute error 



1.9 

3.0 

1.9 


fPercent error = [(calc. — exp.)/exp.j X 100. 

$A11 experimental viscosity values were obtained from Landolt-Bomstein (1955), Lucas (1984), Stephan 
and Lucas (1979). 
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Equation (9-4.16) is similar to an equation proposed by Thodos and coworkers 
(e.g., Yoon and Thodos, 1970). It is interesting to note that, if T r < 1, the f(T r ) in 
brackets in Eq. (9-4.16) is closely approximated by 0.606 T r , that is, 

tj£ = (0.606 T r ) F° P F° q T r < 1 (9-4.20) 

The method of Lucas is illustrated in Example 9-2. 

Example 9-2 Estimate the viscosity of methanol vapor at a temperature of 550 K and 
1 bar by using Lucas’ method. The experimental value is 181 /rP (Stephan and Lucas, 
1979). 

solution From Appendix A, T c = 512.64 K, P c = 80.97 bar, Z c = 0.224, M = 32.042, 
and fi = 1.7 debyes. T r = 550/512.64 = 1.07, and p,,. = 52.46[(1.7) 2 (80.97)/(512.64) 2 ] 

= 4.67 X 10- 2 . From Eq. (9-4.18), 

F° = 1 + (30.55)(0.292 - 0.224) 1 - 72 = 1.30 


With Eq. (9-4.15), 


f = 0.176 


512.64 

(32.042) 2 (80.97) 4 


= 4.70 X 


10-VPL 1 


Then, with Eq. (9-4.16) 


= {(0.807X 1.07)° 618 - 0.357 exp[-(0.449)(1.07)] 
+0.34 exp[-(4.058)(1.07)] + 0.018}(1.30) 

= 0.836 


0.836 

(4.70 x 10 3 ) 


178 n P 


Error = 


178 - 181 
181 


X 100 


-1.7% 


In Table 9-2, experimental viscosities are compared with those computed by Lucas’s 
method. The average absolute error is 3.0%. Even with the correction factor F° p , 
higher errors are noted for polar compounds compared to nonpolar. 

Reichenberg (1971, 1979) has suggested an alternate corresponding states rela¬ 
tion for low-pressure gas viscosity of organic compounds. 

_^(1 + 270 ^) 

' a*[l + (4/r c )][l + 036T r (T r - 1 )] 1/6 T r + 270 

r) is in fjP; M is the molecular weight; T is the temperature; T c is the critical 
temperature, in kelvins; T r is the reduced temperature; and i± r , is the reduced dipole 
moment defined earlier in Eq. (9-4.17). The parameter a* is defined as 

a* = 2n,C, (9-4.22) 

where n, represents the number of groups of the /th type and C, is the group 
contribution shown in Table 9-3. 

The term (1 + 4/7) ) in the denominator of Eq. (9-4.21) may be neglected except 
for treating quantum gases with low values of T c . 


Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com) 
Copyright © 2004 The McGraw-Hill Companies. All rights reserved. 

Any use is subject to the Terms of Use as given at the website. 








VISCOSITY 


VISCOSITY 


TABLE 9-3 Values of the Group Contributions C, for the 
Estimation of a* in Eq. (9-4-22) (Reichenberg, 1971) 

Group 

Contribution C t 

—ch 3 

9.04 

\ 


^CH 2 (nonring) 

6.47 

\ 


CH —(nonring) 

2.67 

\ / 


C (nonring) 

-1.53 

=CH 2 

7.68 

=CH—(nonring) 

5.53 

\ 


C= (nonring) 

1.78 

=CH 

7.41 

=C—(nonring) 

5.24 

\ 

6.91 

^CH,(rmg) 


\ 


CH—(ring) 

1.16 

\ / 


C (ring) 

0.23 

=CH—(ring) 

5.90 

\ 


^C=(rmg) 

3.59 

—F 

4.46 

—Cl 

10.06 

—Br 

12.83 

—OH (alcohols) 

7.96 

\ 


O (nonring) 

3.59 

\ 


C = O (nonring) 

12.02 

—CHO (aldehydes) 

14.02 

—COOH (acids) 

18.65 

—COO—(esters) or HCOO (formates) 

13.41 

—nh 2 

9.71 

\ 


NH (nonring) 

3.68 

/ 



9.13 
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TABLE 9-3 Values of the Group Contributions C, for the 
Estimation of a* in Eq. (9-4-22) (Reichenberg, 1971) 

(i Continued ) 


Group 

Contribution C, 

=N—(ring) 


4.97 

—CN 


18.15 

\ 

S (ring) 


8.86 


A comparison between calculated and experimental low-pressure gas viscosity 
values is given in Table 9-2, and the method is illustrated in Example 9-3. 

Example 9-3 Use Reichenberg’s method to estimate the viscosity of ethyl acetate 
vapor at 125°C and low pressure. The experimental value is reported to be 101 /xP 
(Landolt-Bomstein Tabellen, 1955). 


solution From Appendix A, T c = 523.2 K, M = 88.106, P c = 38.3 bar, and /x = 1.9 
debyes. With Eq. (9-4.17), 


Pr 


(52.46)(1.9) 2 (38.3) 
(523.2) 2 


0.0265 


T r = (125 + 273)/523.2 = 0.761. With Eq. (9-4.22) and Table 9-3 


a* = 2(—CH 3 ) + (—CH 2 ) + (—COO—) 


= (2)(9.04) + 6.47 + 13.41 = 37.96 
With Eq. (9-4.21), 


_ (88.106) 1,2 (125 + 273) 

17 ~ 37.96[ 1 + (0.36)(0.761 X0.761 - l)] 1 ' 6 

(0.761)11 + (270X0.0265 ) 4 ] 

0.761 + (270)(0.0265) 4 ' ^ 

99.4 - 101 

Error = -—- X 100 = -1.5% 


Recommendations for Estimating Low-pressure Viscosities of Pure Gases 

Any of the three estimation methods described in this section may be used with 
the expectation of errors of 0.5 to 1.5% for nonpolar compounds and 2 to 4% for 
polar compounds. Lucas’s method requires as input data T c , P c , and M as well as 
/x and Z c for polar compounds and is easy to apply. At present, it is not suitable 
for highly associated gases like acetic acid, but it could probably be extended by 
multiplication of an appropriate factor as in the Chung, et al. technique. The Chung, 
et al. method requires somewhat more input ( T c , V c , and M and /x, a>, and k for 
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the polar correction). The critical volume is less readily available than the critical 
pressure (See Chap. 2), and the association factor k is an empirical constant that 
must be determined from viscosity data. The method is not suited for quantum 
gases. Reichenberg’s technique requires M, T c , and structural groups as well as yi 
for the polar correction. This method is not suitable for inorganic gases and cannot 
be applied to organic gases for which necessary group contributions have not been 
determined. 


9-5 VISCOSITIES OF GAS MIXTURES 
AT LOW PRESSURES 


The rigorous kinetic theory of Chapman and Enskog can be extended to determine 
the viscosity of low-pressure multicomponent mixtures (Brokaw, 1964, 1965, 1968, 
1965; Chapman and Cowling, 1939; Hirschfelder, et al., 1954; Kestin, et al., 1976). 
The final expressions are quite complicated and are rarely used to estimate mixture 
viscosities. Three simplifications of the rigorous theoretical expressions are de¬ 
scribed below. Reichenberg’s equations are the most complex, but, as shown later, 
the most consistently accurate. Wilke’s method is simpler, and that of Heming and 
Zipperer is even more so. All these methods are essentially interpolative; i.e., the 
viscosity values for the pure components must be available. The methods then lead 
to estimations showing how the mixture viscosity varies with composition. Later 
in this section, two corresponding states methods are described; they do not require 
pure component values as inputs. A compilation of references dealing with gas 
mixture viscosities (low and high pressure) has been prepared by Sutton (1976). 


Method of Reichenberg (1974, 1977, 1979) 

In this technique, Reichenberg has incorporated elements of the kinetic theory ap¬ 
proach of Hirschfelder, et al. (1954) with corresponding states methodology to 
obtain desired parameters. In addition, a polar correction has been included. The 
general, multicomponent mixture viscosity equation is: 


= 2 K, 1 + 2 £ H tj Kj + 2 £ 


j= i 


(9-5.1) 


where r/ m is the mixture viscosity and n is the number of components. With 17 , the 
viscosity of pure i, M, the molecular weight of i, and y, the mole fraction of i in 
the mixture. 


K, = 


yiVi 


(9-5.2) 


y t + v, X y* + (2M,/M,.)] 


Two other component properties used are: 
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U, = 


[1 + 0.3670(7;, - !)]»» F r 


(Tri) 1, 


c, = 


M , 1 


where T ri = 7770 and is a polar correction. 

r ;,: 5 + go / x „-) 7 

7^ 5 [1 + (10 Mr ,.) 7 ] 


(9-5.3) 

(9-5.4) 

(9-5.5) 


Here i± n is the reduced dipole moment of i and is calculated as shown earlier in 
Eq. (9-4.17). For the term H 0 = H jh 


H ,= 


MjMj 


32 (AT, + Mjf 


with 


(C, + C ,) 2 x 


70, = 


[1 + 0.3670,(70, - l )] 1 ' 6 Fj 




(T’ry) 1 ' 


ry (F - T -) 1 

v ci c// 


(9-5.6) 


(9-5.7) 


F RiJ is found from Eq. (9-5.5) with 70 replaced by 70 and /j.,., by = (/x„ /r ) 1/2 . 


For a binary gas mixture of 1 and 2, these equations may be written as: 
Vm = 7^(1 + H\ 2 K\) + K 2 ( 1 + 2 77 12 K\ + 77j 2 F 2 ) 

V, t?, 


Ti + i7i{y 2 77i 2 [3 + (2 M 2 /M 1 )\} 

K = _ME_ 

T 2 + 172 {yj 77 12 [3 + (2 M[/M 2 )]} 


(9-5.8) 

(9-5.9) 


U 1 = 


[1 + 0 . 367 ;,( 7 ;, - i )] 1/6 7 ^ + (io Mfl ) 7 


7,1/2 

J rl 


T^tl + (10 Mrl) 7 ] 


(9-5.10) 


(9-5.11) 


and a comparable expression for U 2 . The meaning of C, and C 2 is clear from Eq. 
(9-5.4). Finally, with 


Trl2 (T cl T c2 ) m 


and LL r i 2 (/-*vi Mi- 2) 1 


(M 1 M 2 /32) 1/2 [1 + 0.36 T rl2 (T rl2 - l )] 1 


Hl2 (M l + M 2 ) 3/2 


x (C, + C 2 ) 2 


M' 2 

T 3 A + (10 /A-12) 7 

1 + (10 /A-12) 7 ] 


(9-5.12) 


To employ Reichenberg’s method, for each component one needs the pure gas 
viscosity at the system temperature as well as the molecular weight, dipole moment, 
critical temperature, and critical pressure. The temperature and composition are state 
variables. 

The method is illustrated in Example 9-4. A comparison of experimental and 
calculated gas-mixture viscosities is shown in Table 9-4. 
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Example 9-4 Use Reichenberg’s method to estimate the viscosity of a nitrogen- 
monochlorodifluoromethane (R-22) mixture at 50°C and atmospheric pressure. The 
mole fraction nitrogen is 0.286. The experimental viscosity is 145 /jlP (Tanaka, et al., 
1977). 


solution The following pure component properties are used: 



n 2 

CHC1F 2 

k 

126.2 

369.28 

P c , bar 

33.98 

49.86 

M, g/mol 

28.014 

86.468 

fj,, debyes 

0 

1.4 

tj, 50°C, juP 

188 

134 


With T = 50°C, T r (N 2 ) = 2.56, and T r (CHC1F 2 ) = 0.875, 


50 + 273.2 
[(126.2)(369.3)] 1/2 


1.497 


tx r (N 2 ) = 0, and from Eq. (9-4.17), 


H r (CHC1F 2 ) = 


(52.46)(1.4) 2 (49.86) 

(369.28) 2 


0.0376 


Since p, rl2 = (/x rl /u. r2 ) 1/2 , then for this mixture, /x rl2 = 0. With Eq. (9-5.11), for CHC1F 2 , 


U (CHC1F 2 ) 


[1 + (0.36)(0.875)(0.875 - 1)] 1/6 (0.875) 3 - 5 + (10) 7 (0.0376) 7 

(0.875) 1/2 X (0.875) 35 [1 + (10) 7 (0.0376) 7 ] 


= 1.063 


and U (N 2 ) = 0.725 


Then, C (N 2 ) = 


(28.014) 1 ' 4 

[(188X0.725)] 1 ' 2 


0.197 


and C (CHC1F 2 ) = 0.256 
Next, 


H(N 2 — CHCIF 2 ) = (0.197 + 0.256) 2 


[(28.014)(86.468)] 1/2 
[32(28.014 + 86.468) 3 ] 1/2 


X 


11 + (0.36)(1.497)( 1.497 - 1)] 1/6 
(1.497) 1/2 


X 1.0 


= 1.237 X 10~ 3 


K (N 2 ) = 


_(0.286)(188)_ 

0.286 + (188)(0.714)(1.237 X 10 3 ){3 + [(2)(86.469)/28.014]} 


29.71 


and K (CHC1F 2 ) = 107.9. Substituting into Eq. (9-5.8), 
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r) m = (29.71)[1 + (1.237 X 10 3 ) 2 ( 107.9) 2 ] + (107.9)[1 

+ (2)( 1.237 X 10 3 )(29.71) + (1.237 X 10 3 ) 2 (29.71) 2 ] 

= 146.2 fiP 
146.2 - 145 

Error = - X 100 = 0.8% 

145 


Method of Wilke (1950) 

In a further simplification of the kinetic theory approach, Wilke (1950) neglected 
second-order effects and proposed: 

V,„ = i (9-5.13) 

' =1 2 y/K 

where 


= 




[8(1 + MilMj)] 112 

<t> j: is found by interchanging subscripts or by 

' " TJ, Mj ' 

For a binary system of 1 and 2, with Eqs. (9-5.13) to (9-5.15), 


(9-5.14) 


(9-5.15) 


yiVi + yiVi 

w + yAn 


where 17 ,,, 
i?i, 172, 

lb, y 2 


viscosity of the mixture 
pure component viscosities 
mole fractions 


and 


Li + ( Vl / V2 y ,2 (M 2 /M 1 y M ] 2 

{ 8[1 + ( Mj / M 2 )]} 1/2 


^*21 


<Al2 


TJlM 2 


(9-5.16) 


Equation (9-5.13), with </_>.. from Eq. (9-5.14), has been extensively tested. Wilke 
(1950) compared values with data on 17 binary systems and reported an average 
deviation of less than 1%; several cases in which 17 ,,, passed through a maximum 
were included. Many other investigators have tested this method (Amdur and Ma¬ 
son. 1958; Bromley and Wilke, 1951; Cheung, 1958; Dahler, 1959; Gandhi and 
Saxena, 1964; Ranz and Brodowsky, 1962; Saxena and Gambhir, 1963, 1963a; 
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Strunk, et al., 1964; Vanderslice, et al. 1962; Wright and Gray, 1962). In most 
cases, only nonpolar mixtures were compared, and very good results obtained. For 
some systems containing hydrogen as one component, less satisfactory agreement 
was noted. In Table 9-4, Wilke’s method predicted mixture viscosities that were 
larger than experimental for the H 2 —N 2 system, but for H 2 —NH 3 , it underestimated 
the viscosities. Gururaja, et al. (1967) found that this method also overpredicted in 
the H 2 —0 2 case but was quite accurate for the H 2 —C0 2 system. Wilke’s approx¬ 
imation has proved reliable even for polar-polar gas mixtures of aliphatic alcohols 
(Reid and Belenyessy, 1960). The principal reservation appears to lie in those cases 
where A7 >> A /. and Vi » Vr 

Example 9-5 Kestin and Yata (1968) report that the viscosity of a mixture of methane 
and n-butane is 93.35 yu.P at 293 K when the mole fraction of n-butane is 0.303. 
Compare this result with the value estimated by Wilke’s method. For pure methane and 
n-butane, these same authors report viscosities of 109.4 and 72.74 fj. P. 


solution Let 1 refer to methane and 2 to n-butane. M x = 16.043 and M, = 58.123. 
[1 + (109.4/72.74) 1/2 (58.123/16.043) 1 ' 4 ] 2 


4* 12 


{8[ 1 + (16.043/58.123)]} 1 


= 2.268 


72.74 16.043 

4> 2X = 2.268 _ = 0.416 


V m = 


109.4 58.123 
(0.697)(109.4) 


(0.303)(72.74) 


0.697 + (0.303)(2.268) 0.303 + (0.697)(0.416) 

= 92.26 yuP 

^ 92.26 - 93.35 

Error = - ^ ^ - X 100 = -1.2% 


Herning and Zipperer (1936) Approximation of <^>, 7 

As an approximate expression for <£,- ; of Eq. (9-5.14) the following is proposed 
(Herning and Zipperer, 1936): 



(9-5.17) 


When Eq. (9-5.17) is used with Eq. (9-5.16) to estimate low-pressure binary gas 
mixture viscosities, quite reasonable predictions are obtained (Table 9-4) except for 
systems such as H 2 — Ni l,. The technique is illustrated in Example 9-6. Note that 
Examples 9-5 and 9-6 treat the same problem; each provides a viscosity estimate 
close to the experimental value. But the (f> 12 and r/) 2l values employed in the two 
cases are quite different. Apparently, multiple sets of and r/^,- work satisfactorily 
in Eq. (9-5.13). 


Example 9-6 Repeat Example 9-5 by using the Herning and Zipperer approximation 
for 
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solution As before, with 1 as methane and 2 as n-butane 

( 58.123 


^12 


= 1.903 


V„, = 


v 16.043) 

<I>21 ~ 4> a = °-525 

(0.697)(109.4) (0.303)(72.74) 


0.697 + (0.303)(1.903) 0.303 + (0.697)(0.525) 


= 92.82 p P 


92.82 - 93.35 

Error = -- X 100 = -0.6% 


Corresponding States Methods 

In this approach, one estimates pseudocritical and other mixture properties (See 
Sec. 5-3) from pure component properties, the composition of the mixture, and 
appropriate combining and mixing rules. 


Lucas (1980, 1983, 1984a) Rules 

Lucas (1980, 1983, 1984a) defined mixture properties as shown below for use in 
Eqs. (9-4.15) through (9-4.19). 


T = y V T 

cm ■ i ci 

i 

(9-5.18) 

2 y,z a 


p _ DT 1 

1 cm cm \ 1 T , 

Z y,K, 

i 

(9-5.19) 

M m = E yMi 

i 

(9-5.20) 

F°Pm = S ViPpi 

i 

(9-5.21) 

FQm ~ ( X yi F °Qi) A 

(9-5.22) 


and, letting the subscript H denote the mixture component of highest molecular 
weight and L the component of lowest molecular weight. 


A = 


0.01 



M 

for —— > 9 and 0.05 < y H < 0.7; 

Mr 


(9-5.23) 


otherwise, A = 1 

The method of Lucas does not necessarily lead to the pure component viscosity 
Tj, when all y, = 0 except y t = 1. Thus the method is not interpolate in the same 
way as are the techniques of Reichenberg, Wilke, and Herning and Zipperer. Nev- 


Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com) 
Copyright © 2004 The McGraw-Hill Companies. All rights reserved. 

Any use is subject to the Terms of Use as given at the website. 







VISCOSITY 


9.24 CHAPTER NINE 

ertheless, as seen in Table 9-4, the method provides reasonable estimates of rj m in 
most test cases. 

Example 9-7 Estimate the viscosity of a binary mixture of ammonia and hydrogen at 
33°C and low pressure by using the Lucas corresponding states method. 


solution Let us illustrate the method for a mixture containing 67.7 mole percent am¬ 
monia. We use the following pure-component values: 



Ammonia 

Hydrogen 

To K 

405.50 

33.2 

P c , bar 

113.5 

13.0 

V c , cm 3 /mol 

72.5 

64.3 

z c 

0.244 

0.306 

M 

17.031 

2.016 

fj., debyes 

1.47 

0 

T r 

0.755 

9.223 


Using Eqs. (9-5.18) to (9-5.20), T cm = 285.2 K, P cm = 89.6 bar, and M m = 12.18. 
From these values and Eq. (9-4.15), = 6.46 X 10 3 (/xP) '. With Eq. 

(9-4.17), tx r (NH 3 ) = 7.825 X 10~ 2 and /x r (H 2 ) = 0. Then, with Eq. (9-4.18), 

F° P (NH 3 ) = 1 + 30.55(0.292 - 0.244) 1 - 72 |0.96 + 0.1(0.755 - 0.7)| = 1.159 

K (H 2 )= 1.0 

F" m = (1.159)(0.677)+ (1)(0.323) = 1.107 

For the quantum correction, with Eq. (9-5.23), since M H /M L = 17.031/2.016 = 
8.4 < 9, then A = 1. F° Q (NH 3 ) = 1.0, and with Eq. (9-4.19), 

F° q (H 2 ) = (1.22)(0.76) o l5 { 1 + 0.00385[(9.209 - 12) 2 ] 1/2016 
X sign (9.209 - 12)} 

= (1.171)[1 + (0.01061)(— 1)] = 1.158 


F" = (1.158X0.323) + (1X0.677) = 1.051 


Next, from Eq. (9-4.16) with T rm = (33 + 273.2)/285.3 = 1.073 
■nJL = (0.645)(1.107)(1.051) = 0.750 
0.750 


V,„ 


6.46 X 10- 3 


= 116.1 /x P 


The experimental value is 120.0 /xP; thus 
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Error = 116 '* * 2Q,Q X 100 = -3.2% 

120.0 

The viscosity of the ammonia-hydrogen mixture at 33°C is line 3 in Fig. 9-3. 


Chung, et al. rules (1984, 1988) 

In this case, Eq. (9-3.9) is employed to estimate the mixture viscosity with, however, 
a factor F cm as used in Eq. (9-4.10) to correct for shape and polarity. 


26.69 F cm {M m TY 


im 2 r\ 

o-fj-k 

where fl„ = In the Chung, et al. approach, the mixing rules are: 


= n.v„ 


T"-/, 


j 1 * = 


M... = 


‘ J 

T 


(e/k) m 

Z X y,,/A)<r, ! , 

j j 

O-m 

'Z 2 y,y/e e /A:)o-?.M !/ 2 


.2 

in''m 


( E /£)„,% 

'Z Z 


(9-5.24) 

(9-5.25) 

(9-5.26) 

(9-5.27) 

(9-5.28) 

(9-5.29) 


Mm °"i 


and the combining rules are: 


* j 

Z Z yjjX, 


y-.yjFF] 


= 0.809E1' 3 


cr, = £„(<r,rr) 


% 

k 

£v 

A: 


\ k k 


1.2593 

1/2 


(9-5.30) 

(9-5.31) 

(9-5.32) 

(9-5.33) 

(9-5.34) 

(9-5.35) 

(9-5.36) 
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", + <Oj 
" y = 2 

(9-5.36) 

Kii = K , 

(9-5.37) 

X,j = (k,Kj )' /2 

(9-5.38) 

2MM, 

M = 1 

,J M, + M 

(9-5.39) 


(9-5.40) 


£ij and 'Cy are binary interaction parameters which are normally set equal to unity. 
The F cm term in Eq. (9-5.24) is defined as in Eq. (9-4.11). 

F cm = 1 - 0.275<u„, + 0.059035 + K m (9-5.41) 

where fj, rm is as in Eq. (9-4.12) 


131-3 /r„, 

(V cm T C J 112 


(9-5.42) 


V cm = (a-,,,/0.809) 3 


(9-5.43) 


T cm = 1.2593 fj] (9-5.44) 

\ IV / m 

In these equations, T c is in kelvins, V c is in cm 3 /mol and /jl is in debyes. 

The rules suggested by Chung, et al. are illustrated for a binary gas mixture in 
Example 9-8. As with the Lucas approach, the technique is not interpolative be¬ 
tween pure component viscosities. Some calculated binary gas mixture viscosities 
are compared with experimental values in Table 9-4. Errors vary, but they are 
usually less than about ±5%. 


Example 9-8 Use the Chung, et al. method to estimate the low-pressure gas viscosity 
of a binary of hydrogen sulfide and ethyl ether containing 20.4 mole percent H,S. The 
temperature is 331 K. 


solution The properties listed below are from Appendix A, the problem statement, 
and Table 9-1: 



Hydrogen sulfide 

Ethyl ether 

T c , K 

373.4 

466.70 

V c , cm 3 /mol 

98 

280 

(X) 

0.090 

0.281 

fj., debyes 

0.9 

1.3 

K 

0 

0 

M, g/mol 

34.082 

74.123 

y 

0.204 

0.796 


From Eqs. (9-5.32) and (9-5.33), 
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ct(H 2 S) = (0.809)(98) 1/3 = 3.730 A 
cr(EE) = 5.293 A 
o-(H 2 S-EE) = 4.443 A 

Then, with Eq. (9-5.25), 

ct 3 = (0.204) 2 (3.730) 3 + (0.796) 2 (5.293) 3 + (2)(0.204)(0.796)(4.443) 3 = 124.58 A 3 
From Eqs. (9-5.34) and (9-5.35) 


« = fH Si 373 ' 4 
k 2 1.2593 

- (EE) = 370.6 K 
k 


= 296.5 K 


7 (H 2 S—EE) = 331.5 K 
k 


Then, with Eq. (9-5.27), 



[(0.204) 2 (296.5)(3.730) 3 + (0.796) 2 (370.6)(5.293) 3 
+ (2)(0.204)(0.796)(331.5)(4.443) 3 ] /124.58 = 360.4 K 


With Eqs. (9-5.28) and (9-5.40), 

M,„ = ({0.204) 2 (296.5)(3.730) 2 (34.082) 1/2 + (0.796) 2 (370.6)(5.293) 2 (74.123) 1/2 
+ (2)(0.204)(0.796)(331.5)(4.443) 2 [(2)(34.082)(74.123)/(34.082 
+ 74.123)] 1/2 } /(360.4)( 124.58) 2 ' 3 ) 2 = 64.44 g/mol 

With Eq. (9-5.29), 

w m = {(0.204) 2 (0.090)(3.730) 3 + (0.796) 2 (0.281)(5.293) 3 

+ (2)(0.204)(0.796)[(0.090 + 0.281)/2](4.443) 3 }/124.58 = 0.256 


and with Eq. (9-5.30), 

l4, = {[(0.204) 2 (0.9) 4 /(3.730) 3 ] + [(0.796) 2 (1.3) 4 /(5.293) 3 ] 

+ [(2)(0.204)(0.796)(0.9) 2 ( 1.3) 2 /(4.443) 3 ]) (124.58) = 2.218 
fi. m = 1.22 debyes 
so, with Eqs. (9-5.42) to (9-5.44), 


V cm 

T cm 


(124.58) 


= 235.3 cm 3 /mol 


(0.809) 3 
= (1.2593X360.4) = 453.9 K 




(131.3)(1.22) 

[(235.3)(453.9)] 1/2 


= 0.490 


Since K m , = 0, with Eq. (9-5.41), 
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F cm = 1 - (0.275)(0.256)+ (0.059035)(0.490) 4 = 0.933 


Using T* from Eq. (9-5.26) [= 331/360.4 = 0.918] and Eq. (9-4.3), 0„ = 1.664. 
Finally, with Eq. (9-5.24), 


Vm 


(26.69)(0.933)[(64.44)(331)] 1/2 
(124.58) 2/3 ( 1.664) 


= 87.6 ,LiP 


The experimental value is 87 /xP (Table 9-4). 

o n fi ot 

Error = ——-X 100 = 0.4% 


Discussion and Recommendations to Estimate the Low-pressure Viscosity 
of Gas Mixtures 

As is obvious from the estimation methods discussed in this section, the viscosity 
of a gas mixture can be a complex function of composition. This is evident from 
Fig. 9-3. There may be a maximum in mixture viscosity in some cases, e.g., system 
3, ammonia-hydrogen. However, cases of a viscosity minimum have been reported. 



0 0.2 0.4 0.6 0.8 1.0 

Mole fraction first component 


FIGURE 9-3 Gas mixture viscosities. 


No. 

System 

Reference 

1 

Hydrogen sulfide-ethyl ether 

Pal and Bhattacharyya (1969) 

2 

Methane-n-butane 

Kestin and Yata (1968) 

3 

Ammonia-hydrogen 

Pal and Barua (1967) 

4 

Ammonia-methyl amine 

Burch and Raw (1967) 

5 

Ethylene-ammonia 

Trautz and Heberling (1931) 
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Behavior similar to that of the ammonia-hydrogen case occurs most often in polar- 
nonpolar mixtures in which the pure component viscosities are not greatly different 
(Hirschfelder, et ah, 1960; Rutherford, et al., 1960). Maxima are more pronounced 
as the molecular weight ratio differs from unity. 

Of the five estimation methods described in this section, three (Heming and 
Zipperer, Wilke, and Reichenberg) use the kinetic theory approach and yield inter- 
polative equations between the pure component viscosities. Reichenberg’s method 
is most consistently accurate, but it is the most complex. To use Reichenberg’s 
procedure, one needs, in addition to temperature and composition, the viscosity, 
critical temperature, critical pressure, molecular weight, and dipole moment of each 
constituent. Wilke’s and Herning and Zipperer’s methods require only the pure 
component viscosities and molecular weights; these latter two yield reasonably 
accurate predictions of the mixture viscosity. 

Arguing that it is rare to have available the pure gas viscosities at the temperature 
of interest, both Lucas and Chung, et al. provide estimation methods to cover the 
entire range of composition. At the end points where only pure components exist, 
their methods reduce to those described earlier in Sec. 9-3. Although the errors 
from these two methods are, on the average, slightly higher than those of the 
interpolative techniques, they are usually less than ±5% as seen from Table 9-4. 
Such errors could be reduced even further if pure component viscosity data were 
available and were employed in a simple linear correction scheme. For example, if 
the pure component viscosity predictions are too high, the mixture prediction would 
be improved if it were lowered by composition-averaged error of the pure com¬ 
ponent predictions. 

An estimation method recently proposed by Davidson (1993) was judged as 
effective as those discussed in this section. 

It is recommended that Reichenberg’s method [Eq. (9-5.8)] be used to calculate 
? 7 m if pure component viscosity values are available. Otherwise, either the Lucas 
method [Eq. (9-4.16)] or the Chung, et al. method LEq. (9-5.24)] can be employed 
if critical properties are available for all components. 


9-6 EFFECT OF PRESSURE ON THE VISCOSITY 
OF PURE GASES 


Figure 9-4 shows the viscosity of carbon dioxide ( T c = 304.1 K and P c = 73.8 
bar) as a function of temperature and pressure. In some ranges ( T r >1.5 and P r < 
2), pressure has little effect on viscosity. But when 1 < T r < 1.5 and when P > 
P c , pressure has a strong effect on viscosity as can be seen by the nearly vertical 
isobars in this region of Fig. 9-4. Figure 9-4 shows isobars as a function of tem¬ 
perature, while Fig. 9-5 shows isotherms as a function of pressure for nitrogen 
(T c = 77.4 K, P c = 33.9 bar). Lucas (1981, 1983) has generalized the viscosity 
phase diagrams (for nonpolar gases) as shown in Fig. 9-6. In this case, the ordinate 
is 17 (; and the temperatures and pressures are reduced values. £ is the inverse reduced 
viscosity defined earlier in Eq. (9-4.15). 

At the critical point, the viscosity diverges so that its value is larger than would 
otherwise be expected. However, this effect is much smaller for viscosity than for 
thermal conductivity (see Fig. 10-5). Whereas the thermal conductivity can increase 
by a factor of two near the critical point, the increase in viscosity is on the order 
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FIGURE 9-4 Viscosity of carbon dioxide. (Stephan and Lucas, 1979) 


of 1%. In fact Vesovic, et al. (1990) state that for carbon dioxide, “the viscosity en¬ 
hancement is less than 1 % at densities and temperatures outside the range bound¬ 
ed approximately by 300K < T < 310K and 300 kg m 3 < p < 600 kg m 3 .” 

In Fig. 9-6, the lower limit of the P r curves would be indicative of the dilute- 
gas state, as described in Sec. 9-4. In such a state, rj increases with temperature. 
At high reduced pressures, we see there is a wide range of temperatures where 17 
decreases with temperature. In this region the viscosity behavior more closely sim¬ 
ulates a liquid state, and, as will be shown in Sec. 9-10, an increase in temperature 
results in a decrease in viscosity. Finally, at very high-reduced temperatures, a 
condition again results in which pressure has little effect and viscosities increase 
with temperature. 

The temperature-pressure region in Fig. 9-4 where viscosity changes rapidly 
with pressure is the very region where density also changes rapidly with pressure. 
Figure 9-7 shows a plot of the residual viscosity as a function of density for n- 
butane. A smooth curve results even though values over a range of temperatures 
are shown. This suggests that density is an important variable when describing 
viscosity behavior at high pressures, and several of the correlations presented in 
this section take advantage of this importance. 
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Enskog Dense-gas Theory 

One of the very few theoretical efforts to predict the effect of pressure on the 
viscosity of gases is due to Enskog and is treated in detail by Chapman and Cowling 
(1939). The theory has also been applied to dense gas diffusion coefficients, bulk 
viscosities, and, for monatomic gases, thermal conductivities. The assumption is 
made that the gas consists of dense, hard spheres and behaves like a low-density 
hard-sphere system except that all events occur at a faster rate due to the higher 
rates of collision (Alder, 1966; Alder and Dymond, 1966) The increase in collision 
rate is proportional to the radial distribution function Mr The Enskog equation for 
shear viscosity is 


— = ^ 1 + 0.8 b oP + 0.761 y(b oP ) 2 (9-6.1) 

V" 

where 17 = viscosity, pP 

rf = low-pressure viscosity, pP 

b a = excluded volume = 2/3 t tN 0 ct 3 , cm 3 /mol 
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FIGURE 9-6 Generalized phase diagram for gas viscosity. (Lucas, 1981, 
1983) 


N a = Avogadro’s number 
cr = hard-sphere diameter, A 
p = molar density, mole/cm 3 

is the radial distribution function at contact and can be related to an equation of 
state by 


'E = Z , 1 (9-6.2) 

(>b„ 

where Z is the compressibility factor. 

Dymond among others (Assael, et al., 1996; Dymond and Assael, 1996) has 
continued efforts to modify the hard sphere approach in order to predict transport 
properties and has shown that viscosities of dense fluids can be correlated by the 
universal equation 
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0 0.10 0.20 0.30 0.40 0.50 0.60 0.70 

Density, g/cm 3 

FIGURE 9-7 Residual //-butane viscosity as a function of density. (Dolan et al., 1963) 


logio Vr = 1-0945 - 


9.26324 71.038 

V, + V? 


301.9012 797.69 

v\ + v? 


1221.977 987.5574 

+ 




Vr 


319.4636 

VI 


(9-6.3) 


where V r = VI V a and V a is a close-packed volume. 17 , is a reduced viscosity defined 
by 


Vr = 6.619 X 10 5 


V V 2 ' 3 

R v (MT) m 


(9-6.4) 


where 17 is in Pa-s, V is in cm 3 /mol, M is g/mol and T is in K. A’ 7( is a parameter 
that accounts for deviations from smooth hard spheres. The two parameters, V a and 
R v are compound specific and are not functions of density. V a is a function of 
temperature as is R for «-alcohols (Assael, et al., 1994). For n-alkanes (Assael, et 
al., 1992a, 1992b), aromatic hydrocarbons (Assael, et al., 1992c), refrigerants (As¬ 
sael, et al., 1995), and a number of other compounds (Assael, et al. 1992; Bleazard 
and Teja, 1996), R y has been found to be independent of temperature. In theory, 
the two parameters V„ and R 7j could be set with two experimental viscosity-density 
data, but in practice, Eq. (9-6.3) has been used only for systems for which extensive 
data are available. It has been applied to densities above the critical density and 
applicability to temperatures down to T r ~ 0.6 has been claimed. Values of V a and 
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R v at 298 K and 350 K (as well as R x which is discussed in Chap. 10) for 16 fluids 
have been calculated with equations given in the above references and are shown 
in Table 9-5. 

Xiang, et al. (1999) have recently extended Eq. (9-6.1) to cover the entire fluid 
range by introducing a crossover function between the low pressure limit and the 
high pressure limit. With a single equation, they fit low pressure viscosities to within 
4% and liquid and high pressure viscosities generally to within 10% for 18 pure 
fluids. Their equation requires the density, critical properties, the acentric factor and 
values for cr and sllc. 


Reichenberg Method (1971, 1975, 1979) 

In this case, the viscosity ratio 17 / 17 " is given by Eq. (9-6.5) 

a pm 

yf ^ BP r + (l + CP?)- 1 


(9-6.5) 


The constants, A, B, C, and D are functions of the reduced temperature T r as shown 
below, and rf is the viscosity of the gas at the same T and low pressure. 


A = ^expa 2 r; B = A((3J r - p 2 ) 


C = — exp y 2 T c r D = — exp S 2 T d r 


TABLE 9-5 Typical values of V„, R v , and R x . 



K, 

cm 3 /mol 





298 K 

350 K 

298 K 

350 K 

298 K 

350 K 

Methane 

17.9 

17.3 

1.00 

1.00 

1.16 

1.16 

n-Butane 

84.1 

84.3 

1.08 

1.08 

1.67 

1.67 

n- Decane 

134.3 

130.1 

1.53 

1.53 

3.13 

3.13 

Cyclohexane 

77.4 

75.7 

0.93 

0.93 

1.35 

1.35 

CC1 4 

65.4 

63.8 

1.07 

1.07 

1.57 

1.57 

R134a 

43.7 

41.7 

1.10 

1.10 

1.61 

1.82 

Ethanol 

33.8 

32.7 

4.62 

2.80 

1.42 

1.42 

rr-Hexanol 

91.9 

84.4 

3.00 

3.00 

1.95 

1.95 

Acetic acid 

41.3 

41.2 

0.76 

0.76 

0.92 

0.98 

Butyl ethanoate 

89.8 

88.3 

1.10 

1.10 

2.20 

2.25 

2-Ethoxyethanol 

70.7 

68.0 

1.33 

1.33 

1.65 

1.78 

1.3-Propanediol 

62.2 

59.1 

1.04 

1.04 

1.16 

1.32 

Diethylene glycol 

78.1 

73.6 

1.83 

1.83 

1.65 

1.94 

Diethanolamine 

87.9 

80.6 

2.38 

2.38 

1.37 

1.72 

Triethylamine 

73.8 

66.4 

2.48 

2.48 

3.03 

3.47 

Dimethyl disulfide 

58.3 

57.0 

1.12 

1.12 

1.79 

1.79 
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a, = 1.9824 X 10 - 3 
p 2 = 1.2760 
8, = 2.9496 


“2 

= 5.2683 

a = 

-0.5767 

7i 

= 0.1319 

72 = 

3.7035 

S 2 

= 2.9190 

d = 

-16.6169 


p 1 = 1.6552 

c = -79.8678 


and Q = (1 — 5.655 yi,), where /jL r is defined in Eq. (9-4.17). For nonpolar materials, 
Q — 1.0. Example 9-9 illustrates the application of Eq. (9-6.5), and, in Table 9-7, 
experimental dense gas viscosities are compared to the viscosities estimated with 
this method. Errors are generally only a few percent; the poor results for ammonia 
at 420 K seem to be an anomaly. 


Example 9-9 Use Reichenberg's method to estimate the viscosity of /i-pentane vapor 
at 500 K and 101 bar. The experimental value is 546 jtiP (Stephan and Lucas, 1979). 


solution Whereas one could estimate the low-pressure viscosity of n-pentane at 500 
K by using the methods described in Sec. 9-4, the experimental value is available (114 
/tiP) (Stephan and Lucas, 1979) and will be used. The dipole moment of n-pentane is 
zero, so Q = 1.0. From Appendix A, T c = 469.7 K and P c = 33.7 bar. Thus T r = 
(500/469.7) = 1.065 and P,. = (101/33.7) = 3.00. From the definitions of A, 5, C, 
and D given under Eq. (9-6.5), A = 0.2999, B = 0.1458, C = 1.271, and D = 7.785. 
With Eq. (9-6.5), 

Jl _ (0.2999)(3.00) 3/2 _ 

■q* (0.1458)(3.00) + [1 + (1.271 )(3.00) 7 - 785 ] 1 

7] = (4.56)(114) = 520 fiP 
520 - 546 

Error = - X 100 = -4.7% 

546 

If one refers back to Fig. 9-6, at T r = 1.065 and P r = 3.00, the viscosity is changing 
rapidly with both temperature and pressure. Thus, an error of only 5% is quite 
remarkable. 


Lucas (1980, 1981, 1983) Method 

In a technique which, in some aspects, is similar to Reichenberg’s, Lucas (1980, 
1981, 1983) recommends the following procedure. For the reduced temperature of 
interest, first calculate a parameter Z, with Eq. (9-6.6). 

Zr = = [0.807T? 618 - 0.357 exp(-0.449 T r ) 

+ 0.340 exp(-4.058 T r ) + 0.018] F° P F° Q (9-6.6) 

where 17 " refers to the low-pressure viscosity. Next calculate Z 2 . If 7’, < 1.0 and P r 
< (P vp /P c ), *en 

Z 2 = 0.600 + 0.760 P" + (6.990 Pf - 0.6)(1 - T r ) (9-6.7) 

with a = 3.262 + 14.98 P 5 r 50S 
P = 1.390 + 5.746 P r 

If (1 < T r < 40) and (0 < P r < 100), then 
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Z 2 = iff 


1 + 


aP e r 


bP f r + (1 + cP d r y\ 


(9-6.8) 


where iff is found from Eq. (9-6.6). The term multiplying this group is identical 
to the pressure correction term in Reichenberg’s method, Eq. (9-6.5), but the values 
of the constants are different. 


a = y exp a 2 TJ 
b = a(b 1 T r — b 2 ) 

^1 rp§ 

c = — exp c 2 T° 


d = y exp d 2 T e r 

e = 1.3088 
f = fi exp f 2 Ti 


1.245 X 10 - 3 

a 2 — 

5.1726 

r = 

-0.3286 

1.6553 

b 2 = 

1.2723 



0.4489 

c 2 = 

3.0578 

8 = 

-37.7332 

1.7368 

d 2 = 

2.2310 

E = 

-7.6351 

0.9425 

Si = 

-0.1853 

i = 

0.4489 


After computing Z, and Z 2 , we define 

T = (9-6.9) 

^1 


and the correction factors F P and F Q , 


F P 


1 + (F° P - 1)T - 3 
F° 

r p 


(9-6.10) 


_ 1 + (F° a - l)[Y~ l - (0.007)(ln 7) 4 ] 

"Q T?o ^ 0.1 

r Q 

where F° P and F° Q are low-pressure polarity and quantum factors determined as 
shown in Eqs. (9-4.18) and (9-4.19). Finally, the dense gas viscosity is calculated 
as 


V = 


Z 2 F p Fq 

£ 


(9-6.12) 


where f is defined in Eq. (9-4.15). At low pressures, Y is essentially unity, and 
F p = 1, F q = 1. Also Z 2 then equals iff so 17 —► if, as expected. 
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The Lucas method is illustrated in Example 9-10, and calculated dense gas 
viscosities are compared with experimental data in Table 9-7. In Example 9-10 and 
Table 9-7, the low-pressure viscosity 17 “ was not obtained from experimental data, 
but was estimated by the Lucas method in Sec. 9-4. Except in a few cases, the 
error was found to be less than 5%. The critical temperature, critical pressure, 
critical compressibility factor, and dipole moment are required, as well as the sys¬ 
tem temperature and pressure. 


Example 9-10 Estimate the viscosity of ammonia gas at 420 K and 300 bar by using 
Lucas’s method. The experimental values of 17 and jf are 571 and 146 /jlP (Stephan 
and Lucas, 1979). 


solution For ammonia we use M = 17.031, Z c = 0.244, T c = 405.50 K, P c = 113.53 
bar, and p. = 1.47 debyes. Thus, T r = (420/405.50) = 1.036 and P r = (300/113.53) 
= 2.643. From Eq. (9-4.15), 


f = (0.176) 


405.50 

(17.031) 3 (113.53) 4 


= 4.95 X lO-^pPL 1 


with Eq. (9-4.17), 


fj. r = (52.46) 


F° a = 1.0 


(1.47) 2 (113.53) 


(405.50) 2 


= 7.827 X 10 ~ 2 


and with Eq. (9-4.18), 

F° p = 1 + 30.55(0.292 - 0.244) L72 |0.96 + (0.1)(1.036 - 0.7)| = 1.164 
From Eq. (9-6.6), Z, = v °£ = 0.7259 

0.7258 


V = 


4.96 x 10- 


= 147 pP 


147 - 146 

Error = - X 100 = 0.7% 

146 


The estimation of the low-pressure viscosity of ammonia agrees very well with the 
experimental value. 

Since T r > 1.0, we use Eq. (9-6.8) to determine Z 2 . The values of the coefficients 
are a = 0.1998, b = 8.834 X 10 2 , c = 0.9764, cl = 9.235, e = 1.3088, and 
f = 0.7808. Then, 


_ (0 1998)(2.643) L3O8S _1 

(8.834 x 10- 2 )(2.643) 0 - 7808 + [1 + (0.9764)(2.643) 9235 ] ‘J ' 

= (4.7761(0.7258) 

= 3.466 

with Eqs. (9-6.9) to (9-6.11), 
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Y 


3.466 

0.7258 


4.775 


Fp 


1 + (1.164 - l)(4.775r 3 
1.164 


0.860 


F q = 1.0 

and, with Eq. (9-6.12), 


= (3.466)(0.860)(1.0) 
V ~ 4.96 x 10- 3 


The Reichenberg and Lucas methods employ temperature and pressure as the state 
variables. In most other-dense gas viscosity correlations, however, the temperature 
and density (or specific volume) are used. In those cases, one must have accurate 
volumetric data or an applicable equation of state to determine the dense gas vis¬ 
cosity. Three different methods are illustrated below. 


Method of Jossi, Stiel, and Thodos (Jossi, et al., 1962; Stiel and Thodos, 
1964) 

In this case, the residual viscosity 17 — 17 " is correlated with fluid density. All 
temperature effects are incorporated in the 17 ° term. To illustrate the behavior of the 
17 — 77 " function, consider Fig. 9-7, which shows tj — 17 " for n-butane graphed as 
a function of density (Dolan, et al., 1963). Note that there does not appear to be 
any specific effect of temperature over the range shown. At the highest density, 0.6 
g/cm 3 , the reduced density pi p c is 2.63. Similar plots for many other substances 
are available, for example. He, air, 0 2 , N 2 , CH 4 (Kestin and Leidenfrost, 1959); 
ammonia (Carmichael, et al., 1963; Shimotake and Thodos, 1963); rare gases (Shi- 
motake and Thodos, 1958); diatomic gases (Brebach and Thodos, 1958); sulfur 
dioxide (Shimotake and Thodos, 1963a); C0 2 (Kennedy and Thodos, 1961; Ve- 
sovic, et al., 1990); steam (Kestin and Moszynski, 1959); and various hydrocarbons 
(Carmichael and Sage, 1963; Eakin and Ellington, 1963; Giddings, 1963; Starling, 
et al., 1960; Starling and Ellington, 1964). Other authors have also shown the 
applicability of a residual viscosity-density correlation (Golubev, 1959; Hanley, et 
al., 1969; Kestin and Moszynski, 1959; Rogers and Brickwedde. 1965; Starling, 
1960, 1962). 

In the lossi, Stiel, and Thodos method, separate residual viscosity expressions 
are given for nonpolar and polar gases, but no quantitative criterion is presented to 
distinguish these classes. 


Nonpolar Gases (Jossi, et al., 1962) 

The basic relation is 

[(T 7 - Tf)£ T + 1 ] 1/4 = 1.0230 + 0.23364 Pr + 0.58533 p 2 r 

- 0.40758 p 3 + 0.093324 p 4 r (9-6.13) 
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where tj = dense gas viscosity, pP 

17 ° = low-pressure gas viscosity, pP 
p r = reduced gas density, pi p c = V c /V 

£ T = the group (TJM i P 4 ) 1/6 , where 7’ is in kelvins and P c is in atmospheres, 
(MP )- 1 

M = molecular weight, g/mol 

This relation is reported by Jossi, et al. to be applicable in the range 0.1 < 


Polar Gases (Stiel and Thodos, 1964) 

The relation to be used depends on the reduced density: 

(tj - t f)£ T = 1.656pJ m p r < 0.1 (9-6.14) 

(77 - vf)Z T = 0.0607(9.045p r + 0.63) 1739 0.1 < p r < 0.9 (9-6.15) 

log{4 - log[(Tj - rf)t T ]} = 0.6439 - 0.1005p r - A 0.9 < p r < 2.6 (9-6.16) 

where A = 0 when 0.9 < p r < 2.2 and 

A = (4.75 X 10- 4 )(p 3 - 10.65) 2 when 2.2 < p r < 2.6 (9-6.17) 

and (77 — rj")£ T = 90.0 and 250 at p r = 2.8 and 3.0, respectively. The notation 
used in Eqs. (9-6.14) to (9-6.17) is defined under Eq. (9-6.13). Note that the pa¬ 
rameter £ T is not the same as £ defined earlier in Eq. (9-4.15). 

An example of the Jossi, et al. method is shown below, and calculated dense 
gas viscosities are compared with experimental values in Table 9-7. 

Example 9-11 Use the Jossi, Stiel, and Thodos method to estimate the viscosity of 
isobutane at 500 K and 100 bar. The experimental viscosity is 261 pP (Stephan and 
Lucas, 1979) and the specific volume is 243.8 cm 3 /mol (Waxman and Gallagher, 1983). 

At low pressure and 500 K, 17 ° = 120 pP. 


solution Since isobutane is nonpolar, Eq. (9-6.13) is used. From Appendix A, T c = 
407.85 K, P c = 36.4 bar = 35.9 atm, V c = 262.7 cm 3 /mol, and M = 58.123. Then 




(407.85) 

(58.123) 3 (35.9) 4 


= 3.282 X 10 2 (pP) 1 


The reduced density = p r = VJV = 262.7/243.8 = 1.078. With Eq. (9-6.13) 
[(tj - 120)(3.282 x 10 2 ) + 1] 1/4 = 1.0230 + (0.23364)( 1.078) 

+ (0.58533)(1.078) 2 
-(0.40758)(1.078) 3 
+ (0.093324)(1.078) 4 


= 1.571 


tj = 275 pP 


275 - 261 

Error = - X 100 = 5.4% 

261 
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Chung, et al. Method (1988) 

In an extension of the Chung, et al. technique to estimate low-pressure gas viscos¬ 
ities, the authors began with Eq. (9-3.9) and employed empirical correction factors 
to account for the fact that the fluid has a high density. Their relations are shown 
below. 


V = 


TJ* 


36.344(M7;) 1/2 

Vf 3 


where 77 = 
M = 
T c = 
^ = 


viscosity, pP 
molecular weight, g/mol 
critical temperature, K 
critical volume, cm 3 /mol 


and 


(T*y 

a 


{^[(Ga 1 + E 6 y]} + 7 ],** 


(9-6.18) 


(9-6.19) 


T * and F c are defined as in Eqs. (9-4.9) and (9-4.11). 11,, is found with Eq. (9-4.3) 
as a function of T*, and, with p in mol/cm 3 , 


y 


PK 

6 


(9-6.20) 


G 1 


1 - 0.5 y 

(i - y ) 3 


£ X {[1 - exp(-£4>')]/>'! + £2^! exp(£ 5 v) + E 2 G l 
E 1 E 4 + E 2 + £3 


(9-6.21) 

(9-6.22) 


77 ** = E 7 y 2 G 2 exp[£ g + E^T*)- 1 + £ 10 (£*r 2 ] (9-6.23) 


and the parameters E l to £ 10 are given in Table 9-6 as linear functions of a> (the 
acentric factor), p * [as defined in Eq. (9-4.12)], and the association factor k (see 
Table 9-1). One might note that, at very low densities, y approaches zero, G , and 
G 2 approach unity, and 77 ** is negligible. At these limiting conditions, combining 


TABLE 9-6 Chung, et al. Coefficients to Calculate E t = a, + b t w + c,-p^ + d t K 


i 

a i 

b t 

C i 

d-i 

1 

6.324 

50.412 

-51.680 

1189.0 

2 

1.210 x 10- 3 

-1.154 x 10~ 3 

-6.257 x 10~ 3 

0.03728 

3 

5.283 

254.209 

-168.48 

3898.0 

4 

6.623 

38.096 

-8.464 

31.42 

5 

19.745 

7.630 

-14.354 

31.53 

6 

-1.900 

-12.537 

4.985 

-18.15 

7 

24.275 

3.450 

-11.291 

69.35 

8 

0.7972 

1.117 

0.01235 

-4.117 

9 

-0.2382 

0.06770 

-0.8163 

4.025 

10 

0.06863 

0.3479 

0.5926 

-0.727 
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Eqs. (9-6.18), (9-6.19) and (9-4.9) leads to Eq. (9-4.10), which then applies for 
estimating 17 ". 

The application of the Chung, et al. method is shown in Example 9-12. Some 
calculated values of 17 are compared with experimental results in Table 9-7. The 
agreement is quite good and errors usually are below 5%. 


Example 9-12 With the Chung, et al. method, estimate the viscosity of ammonia at 
520 K and 600 bar. The experimental value of 17 is 466 71 P (Stephan and Lucas, 1979). 
At this temperature, 17 ” = 182 pP. The specific volume of ammonia at 520 K and 600 
bar is 48.2 cm 3 /mol (Haar and Gallagher, 1978). 


solution We use T c = 405.50 K, = 72.4 cnT/mol, w = 0.256, M = 17.031, and 
71 = 1.47 debyes. Thus T r = 520/405.50 = 1.282 and p = 1/48.2 = 
2.07 X 10 ~ 2 mol/cm 3 . 

With Eq. (9-4.12), 

= (131.3X1.47) = 

^ [(72.4)(405.50)] 1/2 

and with Eq. (9-4.11), 

F c = 1 — (0.2756)(0.256) + (0.059035)(1.13) 4 = 1.026 


T* = (1.2593X1.282) = 1.615 


and with Eq. (9-4.3), = 1.275. Using Eqs. (9-6.20) and (9-6.21), 


(2.075 x 10~ 2 )(72.4) 


= 0.250 and G, = 2.074 


From Table 9-6, the following coefficients were computed: E l = —65.03, £4= —9.287 
X 10~ 3 , £3 = - 204.3, £4 = 2.575, £ 5 = -1.706, £ 6 = 3.018, £ v = 6.749, E s = 
1.103, E g = —1.552, and £ 10 = 1.124. Then, with Eq. (9-6.22), C 2 = 1.472 and, from 
Eq. (9-6.23), 77 ** = 1.101. Finally, using Eqs. (9-6.19) and (9-6.18), 


77* 


(1.615) 1/2 

1.275 


(1.026)[( 1.472) 1 


+ (3.018X0.250)] + 1.101 = 2.567 


V = 


(2.567)(36.344)[(17.031)(405.50)] 1 ' 


(72.4) 2 ' 


= 446 pP 


446 - 466 

Error = -—- X 100 = -4.2% 

466 


TRAPP Method (Huber, 1996) 

The TRAPP (transport property prediction) method is a corresponding states 
method to calculate viscosities and thermal conductivities of pure fluids and 
mixtures. In its original version (Ely, 1981; Ely and Hanely, 1981), it was also used 
to estimate low pressure values of A and 17 and employed methane as a reference 
fluid. In the most recent version presented below for pure fluids and later in Sec. 
9-7 for mixtures, low pressure values are estimated by one of the methods presented 
earlier in the chapter, propane is the reference fluid, and shape factors are no longer 
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functions of density. Other reference fluids could be chosen and in fact, Huber and 
Ely (1992) use R134a as the reference fluid to describe the viscosity behavior of 
refrigerants. The TRAPP method was originally developed only for nonpolar com¬ 
pounds, but there have been efforts to extend the method to polar compounds as 
well (Hwang and Whiting, 1987). 

In the TRAPP method, the residual viscosity of a pure fluid is related to the 
residual viscosity of the reference fluid, propane: 

V ~ rf = F v \ v K - rf"] (9-6.24) 

The reference fluid values are evaluated at T 0 and density p„, not T and p. In Eq. 
(9-6.24), rf is the viscosity at low pressure, rf is the true viscosity of the reference 
fluid, propane, at temperature T a and density p a . rf° is the low pressure value for 
propane at temperature T a . For propane Younglove and Ely (Zaytsev and Aseyev, 
1992) give 


rf-rf’= G, exp [p° l G 2 + p« 5 (p? - DG 3 ] - Gj 


where pf = pj p K , rf 


rf° is in pPa-s, and 


(9-6.25) 


Gi 

G 2 

g 3 

E l 
E 2 
E 3 
E 4 
E 5 
E 6 

Er 

T a , p a , and f' T| are calculated by 

T a = 

Po = 


= exp(£[ + E 2 /T) 

(9-6.26) 

= E 3 + EJT 15 

(9-6.27) 

= £ 5 + ^ ^ 

5 T T 2 

(9-6.28) 

= -14.113294896 


= 968.22940153 


= 13.686545032 


= -12511.628378 


= 0.0168910864 


= 43.527109444 


= 7659.4543472 


T/f 

(9-6.29) 

ph 

(9-6.30) 

( M /YY- 

U 4.094 ) 

(9-6.31) 


where / and h are equivalent substance reducing ratios and are determined as de¬ 
scribed below. 

If vapor pressure and liquid density information are available for the substance 
of interest, and if T < T c , it is recommended that / be obtained from the equation 
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P 


vp 


p 


S 


nm 

P SR (T 0 ) 


(9-6.32) 


where P vp and p s are the vapor pressure and saturated liquid density at temperature, 
T. P R v {T n ) and p SR (T a ) are for the reference fluid, propane. Because the density and 
vapor pressure of the reference fluid are evaluated at T a = T/f, Eq. (9-6.32) must 
be solved iteratively. Once f is found from Eq. (9-6.32), h is determined from 

h = p SR (T 0 )/p s (9-6.33) 


If T > T c , or if vapor pressure and saturated liquid density information are not 
available, h and f can be calculated by 


f = ^ [1 + (w - &> R )(0.05203 - 0.7498 in T r ] (9-6.34) 

h = - 2 $[l~(<o- o/)(0.1436 - 0.2822 in T r ] (9-6.35) 

Pc Z c 


The application of the TRAPP method is shown in Example 9-13. Some calculated 
values of 17 are compared with experimental results in Table 9-7. Huber (1996) 
gives results of additional comparisons and also suggests methods to improve pre¬ 
dictions if some experimental data are available. 


Example 9-13 Repeat Example 9-11 with the TRAPP method. 


solution From Appendix A, for the reference fluid, propane, T c = 369.83 K, V c = 
200 cm 3 /mol, Z c = 0.276 and w = 0.152. For isobutane, T c = 407.85 K, V c = 259 
cm 3 /mol, Z„ = 0.278 and a> = 0.186. With Eqs. (9-6.34) and (9-6.35) followed by 
Eqs. (9-6.29) to (9-6.31) 


f = T 7 TT 7 T [1 + (0.186 - 0.152)(0.05203 - 0.7498 ln( 1.226))] = 1.099 
369.83 

259 X 0 276 

h = onn 3 n 0-70 l 1 “ (°' 186 “ «.152)(0.1436 - 0.2882 ln(1.226))] = 1.282 

ZUU X U.Z/o 

r o = 500/1.099 = 454.9 K, p„ = 1.282/243.8 = 0.005258 mol/cm 3 
(58.124 x 1.099 


h\, = 


44.094 


(1.282r 2/3 .= 1.020 


For the T a and p a above, Eqs. (9-6.25) to (9-6.28) give 17 " — 17 *° = 14.41 pPa s = 
144.1 pP Eq. (9-6.24) gives 


77 = 120 + 1.020 x 144.1 = 267.0 pP 


Error = 


267 - 261 
261 


X 100 = 2.3% 


Other Corresponding States Methods 

In a manner identical in form with that of Chung, et al., Brule and Starling (1984) 
proposed a different set of coefficients for E 1 to E 10 to be used instead of those in 
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Table 9-6. These are shown in Table 9-8. Note that no polarity terms are included 
and the orientation parameter y has replaced the acentric factor w. If values of y 
are not available, the acentric factor may be substituted. 

The Brule and Starling technique was developed to be more applicable for heavy 
hydrocarbons rather than for simple molecules as tested in Table 9-7. 

Okeson and Rowley (1991) have developed a four-parameter corresponding- 
states method for polar compounds at high pressures, but did not test their method 
for mixtures. 


Discussion and Recommendations for Estimating Dense Gas Viscosities 

Six estimation techniques were discussed in this section. Two (Reichenberg and 
Lucas) were developed to use temperature and pressure as the input variables to 
estimate the viscosity. The other four require temperature and density; thus, an 
equation of state would normally be required to obtain the necessary volumetric 
data if not directly available. In systems developed to estimate many types of prop¬ 
erties, it would not be difficult to couple the PVT and viscosity programs to provide 
densities when needed. In fact, the Brule and Starling method (Brule and Starling, 
1984) is predicated on combining thermodynamic and transport analyses to obtain 
the characterization parameters most suitable for both types of estimations. 

Another difference to be recognized among the methods noted in this section is 
that Reichenberg’s, Jossi, et al.’s and the TRAPP methods require a low-pressure 
viscosity at the same temperature. The other techniques bypass this requirement 
and have imbedded into the methods a low-pressure estimation method; i.e., at low 
densities they reduce to techniques as described in Sec. 9-4. If the Lucas, Chung, 
et al., or Brule-Starling method were selected, no special low-pressure estimation 
method would have to be included in a property estimation package. 

With these few remarks, along with the testing in Table 9-7 as well as evaluations 
by authors of the methods, we recommend that either the Lucas or Chung, et al. 
procedure be used to estimate dense (and dilute) gas viscosities of both polar and 
nonpolar compounds. The Brule-Starling method is, however, preferable when com- 


TABLE 9-8 Brule and Starling Coefficients 
to Calculate E t = a, + b t y 


i 

a i 

b, 

i 

17.450 

34.063 

2 

-9.611 X 10- 4 

7.235 X 10- 3 

3 

51.043 

169.46 

4 

-0.6059 

71.174 

5 

21.382 

-2.110 

6 

4.668 

-39.941 

7 

3.762 

56.623 

8 

1.004 

3.140 

9 

-7.774 x 10- 2 

-3.584 

10 

0.3175 

1.1600 


Note: If y values are not available, use co, the 
acentric factor, or, preferably, obtain from multi¬ 
property analysis by using vapor pressure data. 
(Brule and Starling, 1984) 
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plex hydrocarbons are of interest, but even for those materials, the Chung, et al. 
procedure should be used at low reduced temperatures (T, < 0.5). For nonpolar 
compounds, we recommend the TRAPP method as well as the Lucas or Chung, et 
al. methods. 

Except when one is working in temperature and pressure ranges in which vis¬ 
cosities are strong functions of these variables (See Fig. 9-6), errors for the rec¬ 
ommended methods are usually only a few percent. Near the critical point and in 
regions where the fluid density is approaching that of a liquid, higher errors may 
be encountered. 


9-7 VISCOSITY OF GAS MIXTURES AT 
HIGH PRESSURES 


The most convenient method to estimate the viscosity of dense gas mixtures is to 
combine, where possible, techniques given previously in Secs. 9-5 and 9-6. 


Lucas Approach (Lucas 1980, 1981, 1983) 

In the (pure) dense gas viscosity approach suggested by Lucas, Eqs. (9-6.6) to 
(9-6.12) were used. To apply this technique to mixtures, rules must be chosen to 
obtain T c , P c , M, and /u, as functions of composition. For T c , P c , and M of the 
mixture, Eqs. (9-5.18) to (9-5.20) should be used. The polarity (and quantum) 
corrections are introduced by using Eqs. (9-6.10) and (9-6.11), where F° p and F° Q 
refer to mixture values from Eqs. (9-5.21) and (9-5.22). The parameter Y in Eqs. 
(9-6.10) and (9-6.11) must be based on T cm and P cm . F° P and F° Q , for the pure 
components, were defined in Eqs. (9-4.18) and (9-4.19). 


Chung, et al. (1988) Approach 

To use this method for dense gas mixtures, Eqs. (9-6.18) to (9-6.23) are used. The 
parameters T c , V c , oj. M, /jl and k in these equations are given as functions of 
composition in Sec. 9-5. That is, 


Parameter 

Equations to use 

T cm 

(9-5.44), (9-5.27) 

Yen 

(9-5.43), (9-5.25) 

(i) m 

(9-5.29), (9-5.25) 

Yen 

(9-5.32), (9-5.33), (9-5.25) and (9-5.43) 

A,„ 

(9-5.30), (9-5.25) 

K m 

(9-5.31) 


TRAPP Method (Huber, 1996) 

For gas mixtures at high pressure, the viscosity is determined by a combination of 
the techniques introduced for high-pressure gases (Sec. 9-6) with appropriate mix¬ 
ing rules. The viscosity of the mixture is given by: 
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V,„ - V'L = F wn [rf 


rfo] + A^ENSKOG 


(9-7.1) 


The quantity rf — rf° that appears in Eq. (9-7.1) is for the reference fluid 
propane and is evaluated with Eq. (9-6.25) at T a and p a . The following mixing rules 
are used to determine F llm , T a and p„. 



h,„ = ES yjA 

i j 

(9-7.2) 


f m h m = ES yjjfih 

* J 

(9-7.3) 


jS* 

II 

j? 

oo + 

J? 

(9-7.4) 


fij = (fjj) 1 ' 2 

(9-7.5) 

fi and /;, are determined as in Sec. 9-6. /'„ and p 0 arc calculated by equations similar 
to Eqs. (9-6.29) and (9-6.30): 


II 

e--< 0 

(9-7.6) 


Po = pK = hjv 

(9-7.7) 

Finally, 

F vn , = (44.094)- 1/2 (/7„,)- 2 X X yjj ( f ij M ij ) m (h ij 'f f3 

i j 

(9-7.8) 

where 

2 MM, 

M,, ~ 1 

11 M,- + Mj 

(9-7.9) 

The term, Atj enskog , accounts for size differences (Ely, 1981) and is 

calculated by 


A ENSKOG _ ENSKOG _ ENSKOG 

^ '/ '/m 'lx 

(9-7.10) 

where 

v?™ KOG = 2 Pi Y i + y,yj°fafa 

(9-7.11) 


i i j 


p = density in mols/L, cr is in A, and rj° and rf 


are in pP 


48 

25 77 


Y (6-023 X IQ- 4 ) 


= 9.725 X 10“ 


(Tj = 4.771 h]' 3 

+ ,t j 


(9-7.12) 

(9-7.13) 


Because A-i 7 ENSKOG is a correction based on a hard sphere assumption, Eq. (9-3.8) 
is used to calculate rf. The radial distribution function, g ijt is calculated (Tham 
and Gubbins, 1971) by 
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g,i = a - tr l + 


.e„ + ^es 


(i - f ) 2 


(i - f ) 3 


0 ,. = 


OVCT-, t 


y k °i 


2(T ij E 


ytfli 


(9-7.14) 


(9-7.15) 


f = (6-023 x 1(T 4 ) - P 2 y,rr 3 
6 


n = t, 


1 + -f (6.023 X 10- 4 ) P 2 


15 j 1 \ M i + Mj 


M, 


-x 

(T Q 
V°V 


(9-7.16) 


(9-7.17) 


The n values of /3, are obtained by solving the n linear equations of the form 

S Bgpj = Y, (9-7.18) 


where 


«„ = 2 E % 


\Af f + M t 


, 5 o 2 M, 

+ 3 MJ S,:/ 3 ^ 


(9-7.19) 


In Eq. (9-7.19), is the Kronecker delta function, 1 if i = j, and 0 if i * j. The 
quantity rjf NSKOG that appears in Eq. (9-7.10) is for a pure hypothethical fluid with 
the same density as the mixture and is determined with Eq. (9-7.11) with cr v defined 
by 


= ( 2 X ; 

. * j 

M r = V V yjjMfal 


(9-7.20) 

(9-7.21) 


M t j and ir :j are defined in Eqs. (9-7.9) and (9-7.13). Huber (1996) tested the TRAPP 
method on a number of binary hydrocarbon mixtures over a wide range of densities 
and reports an average absolute error of about 5%, although, in some cases, sig¬ 
nificantly larger deviations were found. The method is illustrated in Example 9-14. 


Example 9-14 Use the TRAPP method to estimate the viscosity of a mixture of 
80 mol % methane (1) and 20 mol % n-decane at 377.6 K and 413.7 bar. Lee, et al. 
(1966) report at these conditions, p = 0.4484 g/cm 3 and T/ exp = 126 

jtiPa • s, although this value is considerably higher than values reported by Knapstad, 
et al. (1990) at similar conditions. 


solution From Appendix A 



M 

T c , K 

V c , cm 3 /mol 

Z c 

o 

ch 4 

16.043 

190.56 

98.6 

0.286 

0.011 

Ck>H-22 

142.285 

617.7 

624 

0.256 

0.490 
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Using Eqs. (9-7.2) to (9-7.9) and the procedure illustrated in Example 10-8 leads to 
= 0.9819, h m = 0.8664, T„ = 384.6 K, Po = 9.408 mol/L, and rf - rf° = 51.72 
pPa • s. Equation (9-7.8) gives F vm = 1.260. Calculation of Att ™ 8 * 700 requires the 
application of the method described in Eqs. (9-7.10) through (9-7.21). Intermediate 
results include p = 10.86 mol/L, f = 0.3668 with other values shown below: 


‘j 

A 

V Ip A Pa • s 

Sij 

X 10 4 

11 

3.716 

15.0 

2.696 

210.7 

12 

5.314 

9.86 

3.085 

-6.076 

22 

6.913 

13.0 

3.873 

28.18 


When F, = 2.014 and Y 2 = 0.5627 are used in Eq. (9-7.18), this equation is written 
for each of the two components and solved to give = 102.0 and /3, = 221.6. Finally, 
Eq. (9-7.11) gives t)™ 8KO ° = 911 pP. From Eqs. (9-7.20) and (9-7.21), a x = 4.548 A 
and M v = 47.50. Then, for the hypothetical pure fluid, Eq. (9-7.16) gives 


f = 6.023 X 10- 4 - (10.86)(4.548) 3 = 0.3222 
6 


Eq. (9-7.14), with 8 n = Vi , gives 


1 1 3 X 0.3222 2(0.3222) 2 1 „ 

: — = 2.694 


1 - 0.3222 2 (1 - 0.3222) 2 (1 - 0.3222) 3 4 

Eq. (9-7.17) gives 

Y x = 1 + (6.023 X 10 4 )(10.86) ^ 4 - 548 H 2 - 694 > = 2.339 

From Eq. (9-3.8) 


V° 


^ (47.50 X 377.6 ) 1 ' 2 

26.69 1 -—-— 

(4.548) 2 


172.8 /aP 


For a pure component, Eq. (9-7.19) reduces to B xx = so that B xx = 2.694/172.8 

= 0.01559. Then Eq. (9-7.18) gives j3 x = 2.389/0.01559 = 153.2. Applying Eq. 
(9-7.11) to a pure fluid gives 

tjPnskog = 153 2 x 2.389 + 9.725 X 10“’ (10.86) 2 (4.548) 6 (172.8)(2.694) = 839 ,uP 
With Eq. (9-7.10) 


At/ enskog = 91 i - 839 = 72 pP = 7.2 pPa ■ s. 

Then using Eq. (9-7.1) 

rj m = 10.2 + 1.260 x 51.72 + 7.2 = 82.6 pPa ■ s. 

82.6 - 126 

Error = - X 100 = -34% 

126 

In the above example, for the low pressure contribution of 10.2 pPa-s, Eq. 
(9-3.9) was used for methane, (9-4.16) was used for decane, and (9-5.16) was used 
for the mixture. Although the TRAPP prediction in Example 9-14 was lower than 
the experimental value reported by Lee, et al. (1966), TRAPP predictions are con¬ 
siderably higher than the experimental values reported in (Rnapstad, et al., 1990) 
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for the same system at similar conditions. High-quality data for mixtures with dif¬ 
ferent size molecules at high pressures are limited. This in turn limits the ability 
to evaluate models for this case. 


Discussion 

Both the Lucas and Chung, et al. methods use the relations for the estimation of 
dense gas viscosity and apply a one-fluid approximation to relate the component 
parameters to composition. The TRAPP method uses the term Ai 7 ENSKOG to improve 
the one-fluid approximation. In the Lucas method, the state variables are T, P, and 
composition, whereas in the TRAPP and Chung, et al. procedures, T, p, and com¬ 
position are used. 

The accuracy of the Lucas and Chung, et al. forms is somewhat less than when 
applied to pure, dense gases. Also, as noted at the end of Sec. 9-6, the accuracy is 
often poor when working in the critical region or at densities approaching those of 
a liquid at the same temperature. The TRAPP procedure can be extended into the 
liquid region. The paucity of accurate high-pressure gas mixture viscosity data has 
limited the testing that could be done, but Chung, et al. (1988) report absolute 
average deviations of 8 to 9% for both polar and nonpolar dense gas mixtures. A 
comparable error would be expected from the Lucas form. The TRAPP method 
gives similar deviations for nonpolar mixtures, but has not been tested for polar 
mixtures. Tilly, et al. (1994) recommended a variation of the TRAPP method to 
correlate viscosities of supercritical fluid mixtures in which various solutes were 
dissolved in supercritical carbon dioxide. 

As a final comment to the first half of this chapter, it should be noted that, if 
one were planning a property estimation system for use on a computer, it is rec¬ 
ommended that the Lucas, Chung, et al., or Brule and Starling method be used in 
the dense gas mixture viscosity correlations. Then, at low pressures or for pure 
components, the relations simplify directly to those described in Secs. 9-4 to 9-6. 
In other words, it is not necessary, when using these particular methods, to program 
separate relations for low-pressure pure gases, low-pressure gas mixtures, and high- 
pressure pure gases. One program is sufficient to cover all those cases as well as 
high-pressure gas mixtures. 


9-8 LIQUID VISCOSITY 


Most gas and gas mixture estimation techniques for viscosity are modifications of 
theoretical expressions described briefly in Secs. 9-3 and 9-5. There is no compa¬ 
rable theoretical basis for the estimation of liquid viscosities. Thus, it is particularly 
desirable to determine liquid viscosities from experimental data when such data 
exist. Viswanath and Natarajan (1989) have published a compilation of liquid vis¬ 
cosity data for over 900 compounds and list constants that correlate these data. 
Liquid viscosity data can also be found in Gammon, et al. (1993-1998), Riddick, 
et al. (1986), Stephan and Lucas (1979), Stephen and Hildwein (1987), Stephan 
and Heckenberger (1988), Timmermans (1965), and Vargaftik, et al. (1996). Data 
for aqueous electrolyte solutions may be found in Kestin and Shankland (1981), 
Lobo (1990), and Zaytsev and Aseyev (1992). Tabulations of constants have been 
published in Daubert, et al. (1997), Duhne (1979), van Velzen et al. (1972), Yaws, 
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et al. (1976), and Yaws (1995, 1995a) that allow estimations of liquid viscosities. 
When these constants are derived from experimental data they can be used with 
confidence, but sometimes (Yaws, 1995, 1995a) they are based on estimated vis¬ 
cosities, and in such instances, they should be used only with caution. Liquid phase 
viscosity values can also be found in Dean (1999), Lide (1999), and Perry and 
Green (1997). 

The viscosities of liquids are larger than those of gases at the same temperature. 
As an example, in Fig. 9-8, the viscosities of liquid and vapor benzene are plotted 
as functions of temperature. Near the normal boiling point (353.4 K), the liquid 
viscosity is about 36 times the vapor viscosity, and at lower temperatures, this ratio 
increases even further. Two vapor viscosities are shown in Fig. 9-8. The low- 
pressure gas line would correspond to vapor at about 1 bar. As noted earlier in 
Eq. (9-4.20), below T c , low-pressure gas viscosities vary in a nearly linear manner 
with temperature. The curve noted as saturated vapor reflects the effect of the in¬ 
crease in vapor pressure at higher temperatures. The viscosity of the saturated vapor 
should equal that of the saturated liquid at the critical temperature (for benzene, 
T c = 562.0 K). 



FIGURE 9-8 Viscosities of liquid and vapor benzene. (T b = 353.2 K; 1\. = 562.0 K). 
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Much of the curvature in the liquid viscosity-temperature curve may be elimi¬ 
nated if the logarithm of the viscosity is plotted as a function of reciprocal (abso¬ 
lute) temperature. This change is illustrated in Fig. 9-9 for four saturated liquids: 
ethanol, benzene, n-heptane, and nitrogen. (To allow for variations in the temper¬ 
ature range, the reciprocal of the reduced temperature is employed.) Typically, the 
normal boiling point would be at a value of 77 1 ~ 1.5. For temperatures below the 
normal boiling point (7’, 1 > 1.5), the logarithm of the viscosity varies linearly with 
77*. Above the normal boiling point, this no longer holds. In the nonlinear region, 
several corresponding states estimation methods have been suggested, and they are 
covered in Sec. 9-12. In the linear region, most corresponding states methods have 
not been found to be accurate, and many estimation techniques employ a group 
contribution approach to emphasize the effects of the chemical structure on viscos¬ 
ity. The curves in Fig. 9-9 suggest that, at comparable reduced temperatures, vis¬ 
cosities of polar fluids are higher than those of nonpolar liquids such as hydrocar¬ 
bons, which themselves are larger than those of simple molecules such as nitrogen. 
If one attempts to replot Fig. 9-9 by using a nondimensional viscosity such as rjf 
[see, for example, Eqs. (9-4.13) to (9-4.15)] as a function of T r , the separation 
between curves diminishes, especially at T r > 0.7. However, at lower values of T„ 
there are still significant differences between the example compounds. 

In the use of viscosity in engineering calculations, one is often interested not in 
the dynamic viscosity, but, rather, in the ratio of the dynamic viscosity to the 
density. This quantity, called the kinematic viscosity, would normally be expressed 



1.0 1.5 2.0 


1 /T r 

FIGURE 9-9 Viscosities of various liquids as functions of tempera¬ 
ture. (Stephan and Lucas, 1979) 
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in m 2 /s or in stokes. One stoke (St) is equivalent to 1CL 4 m 2 /s. The kinematic 
viscosity v, decreases with increasing temperature in a manner such that In v is 
nearly linear in temperature for both the saturated liquid and vapor as illustrated in 
Fig. 9-10 for benzene. As with the dynamic viscosity, the kinematic viscosities of 
the saturated vapor and liquid become equal at the critical point. 

The behavior of the kinematic viscosity with temperature has led to several 
correlation schemes to estimate v rather than 17 . However, in most instances. In v 
is related to T l rather than T. If Fig. 9-10 is replotted by using T again there is 
a nearly linear correlation with some curvature near the critical point (as there is 
in Fig. 9-9). 

In summary, pure liquid viscosities at high reduced temperatures are usually 
correlated with some variation of the law of corresponding states (Sec. 9-12). At 
lower temperatures, most methods are empirical and involve a group contribution 
approach (Sec. 9-11). Current liquid mixture correlations are essentially mixing 
rules relating pure component viscosities to composition (Sec. 9-13). Little theory 
has been shown to be applicable to estimating liquid viscosities (Andrade, 1954; 
Brokaw, et al., 1965; Brush, 1962; Gemant, 1941; Hirschfelder, et al., 1954). 



280 320 360 400 440 480 520 560 

Temperature, K 


FIGURE 9-10 Kinematic viscosities of saturated liquid and va¬ 
por benzene (T b = 353.2 K; T c = 562.0 K). 
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9-9 EFFECT OF HIGH PRESSURE ON 
LIQUID VISCOSITY 


9.55 


Increasing the pressure over a liquid results in an increase in viscosity. Lucas (1981) 
has suggested that the change may be estimated from Eq. (9-9.1) 


_ 1 + D(AP r /2.118) A 
Vsl 1 + CwAP r 


(9-9.1) 


where 17 = viscosity of the liquid at pressure P 
t) sl = viscosity of the saturated liquid at P vp 

A P r - (P ~ P^)!P C 

co = acentric factor 

A = 0.9991 - [4.674 x Kr 4 /(1.0523T; 003877 - 1.05 1 3)] 

D = [0.3257/(1.0039 - r 7 - 573 ) 0 - 2906 ] - 0.2086 
C = -0.07921 + 2.16167; - 13.40407; + 44.17067 3 - 
96.12097; - 59.81277® + 15.67197? 


84.82917; + 


In a test with 55 liquids, polar and nonpolar, Lucas found errors in the calculated 
viscosities of less than 10%. To illustrate the predicted values of Eq. (9.9.1), Figs. 
9-11 and 9-12 were prepared. In both, 17 / -r] SL was plotted as a function of A P r for 
various reduced temperatures. In Fig. 9-11, co = 0, and in Fig. 9-12, <0 = 0.2. 
Except at high values of T r , 17 / t) sl is approximately proportional to A P r . The effect 
of pressure is more important at the high reduced temperatures. As the acentric 
factor increases, there is a somewhat smaller effect of pressure. The method is 
illustrated in Example 9-15. 


Example 9-15 Estimate the viscosity of liquid methylcyclohexane at 300 K and 500 
bar. The viscosity of the saturated liquid at 300 K is 0.68 cP, and the vapor pressure 
is less than 1 bar. 


solution From Appendix A, T c = 572.19 K, P c = 34.71 bar, and co = 0.235. Thus 
T r = 300/572.19 = 0.524 and A P r = 500/34.71 = 14.4. (P^ was neglected.) Then 


A = 0.9991 


4.674 X 10 - 4 

(1.0523X0.524) 003877 - 1.0513 


0.9822 


D = 


_0.3257_ 

11.0039 - (0.524) 2 - 573 ] 0 - 2906 


- 0.2086 = 0.1371 


C = -0.07921 + (2.1616)(0.524) - (13.4040)(0.524) 2 + (44.1706)(0.524) 3 
- (84.8291)(0.524) 4 + (96.1209)(0.524) 5 - (59.8127)(0.524) 6 
+ (15.6719)(0.524) 7 = 0.0619 


With Eq. (9-9.1), 


rj _ 1 + (0.137X14.4/2.118)°- 9822 _ 

VsL ~ 1 + (0.235)(14.4)(0.0619) “ ' 

17 = (1.57)(0.68) = 1.07 cP 

The experimental value of 17 at 300 K and 500 bar is 1.09 cP (Titani, 1929). 
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0 2 4 6 8 10 12 14 


P,-Pvp=(P-PvpVP c 

FIGURE 9-11 Effect of pressure on the viscosity of liquids a> = 0. 


1.07 - 1.09 

Error = -—-X 100 = -1.8% 

Whereas the correlation by Lucas would encompass most pressure ranges, at 
pressures over several thousand bar the data of Bridgman suggest that the logarithm 
of the viscosity is proportional to pressure and that the structural complexity of the 
molecule becomes important. Those who are interested in such high-pressure 
regions should consult the original publications of Bridgman (1926) or the work of 
Dymond and Assael (See Sec. 9.6, Assael, et al., 1996, or Dymond and Assael. 
1996). 


9-10 EFFECT OF TEMPERATURE ON 
LIQUID VISCOSITY 


The viscosities of liquids decrease with increasing temperature either under isobaric 
conditions or as saturated liquids. This behavior can be seen in Fig. 9-9, where, for 
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0 2 4 6 8 10 12 14 


P,-P, p ={P-P, p VP c 

FIGURE 9-12 Effect of pressure on the viscosity of liquids; &> = 
0 . 2 . 


example, the viscosity of saturated liquid benzene is graphed as a function of tem¬ 
perature. Also, as noted in Sec. 9-8 and as illustrated in Fig. 9-10, for a temperature 
range from the freezing point to somewhere around the normal boiling temperature, 
it is often a good approximation to assume In r\ L is linear in reciprocal absolute 
temperature; i.e.. 


In m. = A + | (9-10.1) 

This simple form was apparently first proposed by de Guzman (1913) (O’Loane, 
1979), but it is more commonly referred to as the Andrade equation (1930, 1934). 
Variations of Eq. (9-10.1) have been proposed to improve upon its correlation ac¬ 
curacy; many include some function of the liquid molar volume in either the A or 
B parameter (Bingham and Stookey, 1939; Cornelissen and Waterman, 1955; Ev- 
ersteijn, et al., 1960; Girifalco, 1955; Gutman and Simmons, 1952; Innes, 1956; 
Marschalko and Barna, 1957; Medani and Hasan, 1977; Miller, 1963, 1963a; Te- 
lang, 1945; and van Wyk, et al., 1940). Another variation involves the use of a 
third constant to obtain the Vogel equation (1921), 
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FIGURE 9-13 Lewis and Squires liquid 
viscosity-temperature correlation. (Lewis 
and Squires, 1934 as adapted ill Gambill, 
1959) 


In Vl = A + (9-10.2) 

Goletz and Tassios (1977) have used this form (for the kinematic viscosity) and 
report values of A, B, and C for many pure liquids. 

Equation (9-10.1) requires at least two viscosity-temperature datum points to 
determine the two constants. If only one datum point is available, one of the few 
ways to extrapolate this value is to employ the approximate Lewis-Squires chart 
(1934), which is based on the empirical fact that the sensitivity of viscosity to 
temperature variations appears to depend primarily upon the value of the viscosity. 
This chart, shown in Fig. 9-13, can be used by locating the known value of viscosity 
on the ordinate and then extending the abscissa by the required number of degrees 
to find the new viscosity. Figure 9-13 can be expressed in an equation form as 

v, 0266 ' = v* 0 - 2661 + (9-10.3) 

where i) L = liquid viscosity at T, cP 

iq K = known value of liquid viscosity at T K , cP 

T and T K may be expressed in either °C or K. Thus, given a value of rp at T K , one 
can estimate values of 17 , at other temperatures. Equation (9-10.3) or Fig. 9-13 is 
only approximate, and errors of 5 to 15% (or greater) may be expected. This method 
should not be used if the temperature is much above the normal boiling point. 

Example 9-16 The viscosity of acetone at 30°C is 0.292 cP; estimate the viscosities 
at -90°C, -60°C, 0°C, and 60°C. 
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solution At -90°C, with Eq. (9-10.3), 


^-0.2661 = ( o.292 )-°- 2661 


-90 - 30 
233 


Vl = 1-7 cP 


For the other cases, 


u 

0 

Vl, cP 

Eq. (9-10.3) 

Vl, cP 

Experimental 

Precent 

error 

-90 

1.7 

2.1 

-19 

-60 

0.99 

0.98 

1 

0 

0.42 

0.39 

8 

60 

0.21 

0.23 

-9 


In summary, from the freezing point to near the normal boiling point, Eq. (9-10.1) 
is a satisfactory temperature-liquid viscosity function. Two datum points are required. 
If only one datum point is known, a rough approximation of the viscosity at other 
temperatures can be obtained from Eq. (9-10.3) or Fig. 9-14. 

Liquid viscosities above the normal boiling point are treated in Sec. 9-12. 


9-11 ESTIMATION OF LOW-TEMPERATURE 
LIQUID VISCOSITY 


Estimation methods for low-temperature liquid viscosity often employ structural- 
sensitive parameters which are valid only for certain homologous series or are found 
from group contributions. These methods usually use some variation of Eq. (9-10.1) 
and are limited to reduced temperatures less than about 0.75. We present two such 
methods in this section. We also describe a technique that employs corresponding 
states concepts. None of the three methods considered is particularly reliable. 


Orrick and Erbar (1974) Method 

This method employs a group contribution technique to estimate A and B in Eq. 
(9-11.1). 


In = A + - 
Pl M T 


(9-11.1) 


where 17 , = liquid viscosity, cP 

p L = liquid density at 20°C, g/cm 3 
M = molecular weight 
T = temperature, K 

The group contributions for obtaining A and B are given in Table 9-9. For liquids 
that have a normal boiling point below 20°C, use the value of p L at 20°C; for liquids 
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TABLE 9-9 Orrick and Erbar (1974) 

Group Contributions for A and B in 

Eq. (9-11.1) 

Group 

A 

B 

Carbon atoms t 

-(6.95 + 0.21A9 

275 + 99 N 

X 

R_ r R 

-0.15 

35 

R 

| 



R—C—R 

X 

-1.20 

400 

K 

Double bond 

0.24 

-90 

Five-membered ring 

0.10 

32 

Six-membered ring 

-0.45 

250 

Aromatic ring 

0 

20 

Ortho substitution 

-0.12 

100 

Meta substitution 

0.05 

-34 

Para substitution 

-0.01 

-5 

Chlorine 

-0.61 

220 

Bromine 

-1.25 

365 

Iodine 

-1.75 

400 

—OH 

-3.00 

1600 

—COO— 

-1.00 

420 

—O— 

-0.38 

140 

1 

—c=o 

-0.50 

350 

—COOH 

-0.90 

770 


f N = number of carbon atoms not including those in other groups shown above. 


whose freezing point is above 20°C, p L at the melting point should be employed. 
Compounds containing nitrogen or sulfur cannot be treated. Orrick and Erbar tested 
this method for 188 organic liquids. The errors varied widely, but they reported an 
average deviation of 15%. This is close to the average value of 16% shown in Table 
9-11 for a more limited test. Since p L in Eq. (9-11.1) is at 20°C and not T, the 
temperature of the liquid, Eq. (9-11.1) is the same form as the Andrade equation, 
Eq. (9-10.1). 


Example 9-17 Estimate the viscosity of liquid n-butyl alcohol at 120°C with the 
Orrick-Erbar method. The experimental value is 0.394 cP. 

solution From Table 9-9 

A = -6.95 - (0.21)(4) - 3.00 = -10.79 
B = 275 + (99)(4) + 1600 = 2271 

From Vargaftik, et al. (1996), at 20°C, p L = 0.8096 g/cm 3 and M = 74.123. Then, 
with Eq. (9-11.1), 
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In 


Vl 

(0.8096X74.123) 


-10.79 


2271 

T 


At T= 120°C = 393 K, r) L = 0.400 cP 


Error = 


0.400 - 0.394 
0.394 


X 100 = 1.5% 


Sastri-Rao Method (1992) 

In this method, the pure liquid viscosity is calculated with the equation 

V = r h ,P^ (9-11.2) 

P vp is the vapor pressure in atmospheres and r/ s is the viscosity at the normal boiling 
point, T b in mPa-s. Below T b , Sastri and Rao determine P with the equation 

In P vp = (4.5398 + 1.0309 In T b ) 

j (9 _ 9 TIT \0 19 \ 

X I 1 - -- — b) -0.38(3 - 2TIT h ) 019 ln(T/T b )J (9-11.3) 

Equation (9-11.3) should be used only when T < T b . Equation (9-11.3) is not 
necessarily the most accurate equation for vapor pressure predictions but should be 
used with Eq. (9-11.2) because the group contributions used to estimate rj,, and N 
have been determined when P vp was calculated with Eq. (9-11.3). r] B is determined 
with the equation 


Vb = 2 A 17 B + 2 A 17 Bear (9-11.4) 

N is determined from 

N = 0.2 + 2 AM + 2 A/V„„ (9-11.5) 

Values for group contributions to determine the summations in Eqs. (9-11.4) and 
(9-11.5) are given in Table 9-10. The contributions of the functional groups to r) B 
and N are generaly cumulative. However, if the compound contains more than one 
identical functional group, its contributions for N should be taken only once unless 
otherwise mentioned. Thus for branched hydrocarbons with multiple >CH- groups, 
N is 0.25. In Table 9-10, the term alicyclic means cycloparaffins and cycloolefins 
and excludes aromatics and heterocyclics. In the contributions of halogen groups, 
“others” means aromatics, alicyclics, and heterocyclics while the carbon groups 
listed are meant for aliphatic compounds. Also for halogens, the values of A N for 
aliphatic, alicyclics and aromatics are not used if other non-hydrocarbon groups are 
present in the cyclic compound (See footnote b in the halogen section of Table 
9-10). For example, the corrections for halogenated pyridines and anilines are given 
in footnote b and are not to be used in conjunction with the corrections listed under 
“aliphatic, alicyclics and aromatics.” Calculation of r\ B and N is illustrated in Ex¬ 
ample 9-18 and typical deviations are shown in Table 9-11. 

Example 9-18 Determine the values of tj b and N to be used in Eq. (9-11.2) for o- 
xylene, ethanol, ethylbenzene, 2-3-dimethylbutane, and o-chlorophenol. 
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TABLE 9-11 Comparison of Calculated and Experimental Viscosities of Liquids 


Percent error* in liquid viscosity 
calculated by the method of 


Compound 

T, K 

V (exp.) 

C p** 

Orrick 

and 

Erbar 

Sastri 

and 

Rao 

Przezdziecki 
and Sridhar 

Acetone 

183 

2.075 

-25 

-3.5 

-11 


213 

0.982 

-6.7 

1.9 

-4.6 


273 

0.389 

-8.3 

3.6 

-2.3 


303 

0.292 

-9.4 

1.6 

-1.2 


333 

0.226 

-8.3 

1.9 

0.2 

Acetic acid 

283 

1.450 

-22 

-15 

8.6 


313 

0.901 

-15 

-15 

0 


353 

0.561 

-9.5 

-17 

-1.3 


383 

0.416 

-5.3 

-16 

0.3 

Aniline 

263 

13.4 

— 

-24 

— 


293 

4.38 

— 

-4.8 

— 


333 

1.520 

— 

8.1 

-49 


393 

0.658 

— 

-9.6 

-33 

Benzene 

278 

0.826 

-45 

-8.5 

1.1 


313 

0.492 

-35 

-6.6 

7.3 


353 

0.318 

-26 

-5.5 

12 


393 

0.219 

-46 

-5.7 

18 


433 

0.156 

-7.1 

-6.4 

23 


463 

0.121 

5.1 

-9.5 

28 

n-Butane 

183 

0.630 

-14 

1.6 

-9 


213 

0.403 

-20 

-2.8 

-8.9 


273 

0.210 

-23 

-0.4 

-5.8 

1-Butene 

163 

0.79 

-22 

0.9 

-13 


193 

0.45 

-20 

-2.9 

-9.6 


233 

0.26 

-18 

-2.5 

-3.3 

n-Butyl alcohol 

273 

5.14 

-2.1 

0.3 

— 


313 

1.77 

-1.6 

-3.4 

— 


353 

0.762 

0.5 

-2 

— 


393 

0.394 

-1.4 

1.4 

— 

Carbon tetrachloride 

273 

1.369 

20 

-4.4 

-24 


303 

0.856 

22 

-2 

-15 


343 

0.534 

20 

-0.1 

-6.7 


373 

0.404 

19 

-0.7 

-2.8 

Chlorobenzene 

273 

1.054 

1.4 

2.7 

-8.3 


313 

0.639 

-0.6 

0.8 

-7 


353 

0.441 

-0.9 

-1.2 

-5.2 


393 

0.326 

-5.1 

-0.9 

-3.8 

Chloroform 

273 

0.700 

40 

7.4 

-11 


303 

0.502 

34 

5.7 

-8.1 


333 

0.390 

27 

3.6 

-7.9 

Cyclohexane 

278 

1.300 

-51 

-29.7 

-38 


333 

0.528 

-38 

-16.4 

-22 

Cyclopentane 

293 

0.439 

-32 

-5.1 

-33 


323 

0.323 

-28 

-7.8 

-29 
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TABLE 9-11 Comparison of Calculated and Experimental Viscosities of Liquids 
(.Continued) 


Percent error* in liquid viscosity 
calculated by the method of 


Compound 

T, K 

V (exp.) 

qP** 

Orrick 

and 

Erbar 

Sastri 

and 

Rao 

Przezdziecki 
and Sridhar 

2,2-Dimethylpropane 

258 

0.431 

-3.5 

-24.3 

20 


283 

0.281 

-0.8 

-15.1 

30 

Ethane 

98 

0.985 

30 

53.6 

-24 


153 

0.257 

-12 

26.4 

-14 


188 

0.162 

-22 

21.3 

-13 

Ethylene chloride 

273 

1.123 

-43 

-20.7 

— 


313 

0.644 

-35 

-15.3 

— 


353 

0.417 

-27 

-8.8 

— 

Ethyl alcohol 

273 

1.770 

27 

-14.1 

— 


313 

0.826 

3.5 

-6 

— 


348 

0.465 

-5.4 

7.7 

— 

Ethyl acetate 

293 

0.458 

-4.2 

-5.5 

-16 


353 

0.246 

0.4 

-1.2 

-5.3 


413 

0.153 

7.4 

-9.7 

-1.8 


463 

0.0998 

27 

-2.2 

4.8 

Ethylbenzene 

253 

1.240 

-2.9 

19.7 

-33 


313 

0.535 

-1.2 

7.3 

-23 


373 

0.308 

-1.7 

0.1 

-16 


413 

0.231 

-1.2 

-1.9 

-13 

Ethyl bromide 

293 

0.395 

27 

0 

-23 


333 

0.269 

32 

2.2 

-17 


373 

0.199 

36 

4.6 

-16 

Ethylene 

103 

0.70 

-25 

-0.7 

25 


133 

0.31 

-27 

-2.3 

-17 


173 

0.15 

-22 

4.8 

-6.4 

Ethyl ether 

273 

0.289 

0 

2.7 

0 


293 

0.236 

0 

1.6 

2.2 


333 

0.167 

2 

-3.5 

4 


373 

0.118 

1 

5.5 

7.4 

Ethyl formate 

273 

0.507 

-18 

6.8 

-16 


303 

0.362 

-17 

6.6 

-11 


328 

0.288 

-16 

7 

-9.6 

n-Heptane 

183 

3.77 

-21 

-33.2 

-1.7 


233 

0.965 

-0.5 

-11.3 

-27 


293 

0.418 

-1.9 

-7.2 

-21 


373 

0.209 

-3.3 

-1 

-17 

n-Hexacontane 

408 

7.305 

18 

37 

— 

(050^122) 

466 

3.379 

-63 

19 

— 

n-Hexane 

213 

0.888 

2.9 

-2.3 

-8.3 


273 

0.381 

-2.4 

-2.4 

-8.2 


343 

0.205 

-4.9 

1.2 

-7.1 
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TABLE 9-11 Comparison of Calculated and Experimental Viscosities of Liquids 
( Continued ) 


Compound 

r, k 

V (exp.) 

qP** 

Percent error* in liquid viscosity 
calculated by the method of 

Orrick Sastri 

and and Przezdziecki 

Erbar Rao and Sridhar 

Isobutane 

193 

0.628 

-23 

1.1 

-37 


233 

0.343 

-25 

-12.2 

-29 


263 

0.239 

-24 

-17.1 

-23 

Isopropyl alcohol 

283 

3.319 

-24 

-16.2 

— 


303 

1.811 

-15 

9.8 

— 


323 

1.062 

-10 

4.5 

— 

Methane 

88 

0.226 

60 

1.5 

-11 


113 

0.115 

23 

-20.1 

-4.3 

2-Methylbutane 

223 

0.550 

-13 

-10.8 

-30 


253 

0.353 

-12 

0.3 

-21 


303 

0.205 

-10 

-1.9 

-12 

n-Pentane 

153 

2.35 

-1 

-1.7 

11 


193 

0.791 

3.8 

0.1 

-7 


233 

0.428 

-3.3 

7.3 

-6 


273 

0.279 

-8.2 

2.9 

-4.7 


303 

0.216 

-11 

0.6 

-4.9 

Phenol 

323 

3.020 

0 

1.5 

-50 


373 

0.783 

37 

12.5 

-5.4 

Propane 

133 

0.984 

-1.5 

45.4 

-23 


193 

0.327 

-22 

14 

-19 


233 

0.205 

-25 

-1.7 

-16 

n-Propyl alcohol 

283 

2.897 

-9.1 

0.8 

— 


313 

1.400 

-9.8 

-6.1 

— 


373 

0.443 

-6.5 

-7.4 

— 

Toluene 

253 

1.070 

-19 

0.4 

-33 


293 

0.587 

-13 

0.7 

-24 


333 

0.380 

-10 

-2.1 

-16 


383 

0.249 

-6.8 

-5.1 

-10 

o-Xylene 

273 

1.108 

3.1 

-5.1 

-5.5 


313 

0.625 

5 

-4.5 

-4.8 


373 

0.345 

3.7 

4.5 

-0.3 


413 

0.254 

3.6 

1.8 

-2.7 

1.9 

m-Xylene 

273 

0.808 

1.1 

-2.6 

1.9 


313 

0.492 

1.4 

-0.4 

1.8 


353 

0.340 

0.3 

2 

2.9 


413 

0.218 

1.4 

3.4 

4.6 


*[(calc — expl/exp] X 100. 

**Data from Aasen, et al. (1990), Amdur and Mason (1958), and Landolt-Bomstein (1955). 


Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com) 
Copyright © 2004 The McGraw-Hill Companies. All rights reserved. 

Any use is subject to the Terms of Use as given at the website. 






VISCOSITY 


9.72 


CHAPTER NINE 


solution o-xylene has 4=CH— (ring, not alicyclic), 2=C< (ring, not alicyclic) and 
2 —CH 3 groups. There is a correction to i 7 Bcor of 0.07 for multiple substitution. With 
values from Table 9-10: 

i|, = 4 x 0.05 - 2 x 0.12 + 2 X 0.105 + 0.07 = 0.24 mPa-s 
N = 0.2 + 0.05 = 0.25 

ethanol has one —CH 3 , one >CH 2 (non-ring) and one —OH. With values from Table 
9-10: 


i) B = 0.105 + 0.615 - 0.092 X 2 + 0.004 x 4 - l0 -°- 58x2 = 0.483 mPa-s 
N = 0.2 + 0.15 = 0.35 

ethylbenzene has 5 =CH— (ring, not alicyclic), one =C< (ring not alicyclic), one 
—CH 3 , and one —CH 2 — (non-ring). There is a branching correction to A N of 0.025. 
With values from Table 9-10: 

= 5 x 0.05 - 0.12 + 0.105 = 0.235 mPa-s 

N = 0.2 + 0.025 = 0.225 

2-3 dimethylbutane has 4 —CH 3 and 2 >CH— (non-ring). The value of AN of 0.05 
is applied only once. With values from Table 9-10: 

t) s = 4 X 0.105 - 2 X 0.11 = 0.2 mPa-s 

N = 0.2 + 0.05 = 0.25 

o-chlorophenol has 4=CH— (ring, not alicyclic), 2=C< (ring, not alicyclic), one — 
Cl attached to an “other”, and one —OH (phenolic). Note that the —Cl contribution 
to AN of 0.025 is not used. Footnote b in the halogen section of Table 9-10 applies 
because of the presence of the non-hydrocarbon —OH group. With values from Table 
9-10: 


r, B = 4 x 0.05 - 2 x 0.12 + 0.17 + 0.27 = 0.4 mPa-s 
N = 0.2 + 0.2 - 0.075 = 0.325 


Przezdziecki and Sridhar (1985) Method 


In this technique, the authors propose using the Hildebrand-modified Batschinski 
equation (Batschinski, 1913; Hildebrand, 1971; Vogel and Weiss, 1981) 


V l = 


Vo 

E(V - V a ) 


(9-11.6) 


where rj, = liquid viscosity, cP 

V = liquid molar volume, cm 3 /mol 

and the parameters E and V D are defined below. 


E = -1.12 + 


12.94 + 0.10 M - 0.23 P c + 0.0424 T fP - 11.58 (T fP /T c ) 


(9-11.7) 
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V a = 0.0085w7; - 2.02 + 


Vn, 

0.342(7^/7;) + 0.894 


where T c 
Pc 
Vc 
M 
Pfp 

(O 

V 


critical temperature, K 
critical pressure, bar 
critical volume, cm 3 /mol 
molecular weight, g/mol 
freezing point, K 
acentric factor 

liquid molar volume at T fp , cm 3 /mol 


(9-11.8) 


Thus, to use Eq. (9-11.6), one must have values for T c , P c , V c , T fP , co, and V m in 
addition to the liquid molar volume V at the temperature of interest. The authors 
recommend that V m and V be estimated from T fp and T by the Gunn-Yamada (1971) 
method. In the Gunn-Yamada method, one accurate value of V is required in the 
temperature range of applicability of Eq. (9-11.6). We define this datum point as 
V R at 7' K : then at any other temperature T, 


V(T) 


f(T) 

HT*) 


V R 


(9-11.9) 


where 


f(T) = //,(! - ojH 2 ) (9-11.10) 

H 1 = 0.33593 - 0.33953 T r + 1.51941 T 2 r - 2.02512 T\ (9-11.11) 
+ 1.114227? 

H 2 = 0.29607 - 0.090457; - 0.048427? (9-11.12) 

Equation (9-11.6) was employed with Eqs. (9-11.7) to (9-11.12) to estimate 
liquid viscosities for the compounds in Table 9-11. The values of T c , P c , V c , T fp , 
and (o , were obtained from Appendix A. The reference volume for each compound 
was calculated from the liquid density datum value given in Appendix A. Large 
errors were noted for alcohols, and those results are not included in the table. For 
other compounds, the errors varied widely and, except for a few materials, the 
technique underestimated the liquid viscosity. Larger errors were normally noted at 
low temperatures, but that might have been expected from the form of Eq. (9-11.6). 
That is, because V„ is of the order of the volume at the freezing point and rj, oc 
(V — V,,) -1 , the estimated value of 17 , becomes exceedingly sensitive to the choice 
of V. This problem was emphasized by Luckas and Lucas (1986), who suggest that 
Eq. (9-11.6) should not be used below T r values of about 0.55. 


Example 9-19 Use the Przezdziecki and Sridhar correlation to estimate the liquid 
viscosity of toluene at 383 K. The experimental value is 0.249 cP (Vargaftik, et af, 
1996). 
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solution From Appendix A, for toluene (slightly different values were used to cal¬ 
culate the results shown in Table 9-11) 


T c = 591.75 K 
P c = 41.08 bar 
V c = 316 cm 3 /mol 
T fp = 178 K 
M = 92.14 g/mol 
w = 0.264 

V L = 106.87 cm 3 /mol at 298.15 K 
With T R = 298.15 K, and with Eqs. (9-11.9) to (9-11.12), 


T * 298.15 
' 591.75 


0.504 


H^T?) = 0.33593 - (0.33953)(0.504) + (1.51941)(0.504) 2 - (2.02512)(0.504) 3 
+ (1.11422)(0.504) 4 = 0.363 

H 2 (T R ) = 0.29607 - (0.09045X0.504) - (0.04842)(0.504) 2 = 0.238 
f(T R ) = 0.363[1 - (0.264X0.238)] = 0.340 
Similarly, 



T, K 

T r 

H\ 

h 2 

f(T) 

T f P 

178 

0.301 

0.325 

0.264 

0.303 

T 

383 

0.647 

0.399 

0.217 

0.376 


0.303 

V= —— (106.87) = 95.2 cm 3 /mol 
Then 0340 

n 77 

V = - -(106.87) = 118.2 cm 3 /mol 

0.340 


This value for V agrees with that given in Vargaftik, et al. (1996). With Eqs. (9-11.7) and 
(9-11.8) 


E = -1.12 + 316/[ 12.94 + (0.10)(92.14) - (0.23)(41.08) + (0.0424)(178) 
- (11.58X178/591.8)] = 17.72 


V a = (0.0085)(0.264)(591.75) - 2.02 + 


_9X2_ 

[(0.342X178/591.8) + 0.894] 


= 94.8 cm 3 /mol 


Then, with Eq. (9-11.6) 


Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com) 
Copyright © 2004 The McGraw-Hill Companies. All rights reserved. 

Any use is subject to the Terms of Use as given at the website. 









VISCOSITY 


VISCOSITY 


9.75 


Vl 


_9T8_ 

17.72(118.2 - 94.8) 


= 0.229 cP 


Error 


0.229 - 0.249 
0.249 


X 100 = -8% 


Other Correlations 

Other viscosity-correlating methods have been proposed, and a number of these are 
summarized in Mehrotra, et al. (1996) and Monnery, et al. (1995). Other recent 
correlations are given in Mehrotra (1991), and the earlier literature was reviewed 
in the 4 lh Edition of this book. 


Recommendations for Estimating Low-temperature Liquid Viscosities 

Three estimation methods have been discussed. In Table 9-11, calculated liquid 
viscosities are compared with experimental values for 36 different liquids (usually 
of simple structure). Large errors may result, as illustrated for all methods. The 
method of Przezdziecki and Sridhar should not be used for alcohols. 

The method of Sastri and Rao assumes that the temperature dependence of 17 , 
is related to the temperature dependence of the vapor pressure, whereas the Orrick 
and Erbar method is slightly modified to include the liquid density. Neither is 
reliable for highly branched structures or for inorganic liquids and the Orrick-Erbar 
method cannot be used for sulfur compounds. Both are limited to a temperature 
range from somewhat above the freezing point to about T r ~ 0.75. Przezdziecki 
and Sridhar’s method employs the Hildebrand equation, which necessitates knowl¬ 
edge of liquid volumes. 

It is recommended that, in general, the method of Sastri and Rao be used to 
estimate low-temperature liquid viscosities. Errors vary widely, but should be less 
than 10 to 15% in most instances. 


9-12 ESTIMATION OF LIQUID VISCOSITY AT 
HIGH TEMPERATURES 


Low-temperature viscosity correlations as covered in Sec. 9-10 usually assume that 
In % is a linear function of reciprocal absolute temperature. Above a reduced 
temperature of about 0.7, this relation is no longer valid, as illustrated in Fig. 
9-10. In the region from about T r = 0.7 to near the critical point, many estimation 
methods are of a corresponding states type that resemble or are identical with those 
used in the first sections of this chapter to treat gases. For this temperature range, 
Sastri (1998) recommends 


In 17 


In y B 

ln(cYi7 B ) 


ln(ai7 s ) 


(9-12.1) 


where 17 is in mPa • s 
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r) B = viscosity at T b in mPa • s, from contributions in Table 9-10 
a = 0.1175 for alcohols and 0.248 for other compounds 

(9-12.2) 

1 1 br 


where T r = TIT C 
Tbr — T b IT c 

Sastri reports average deviations of 10% for T r > 0.9 and 6% for T hr < T r < 0.9. 

Example 9-20 Estimate the saturated liquid viscosity of «-propanol at 433.2 K by 
using Eq. (9-12.1). The experimental value is 0.188 cP. 


solution From Appendix A, T b = 370.93 K and T c = 536.78 K. With contributions 
from Table 9-10, t, b = 0.105 + 0.615 - 0.092 X 3 + 0.004 X 9 - 10- 3XO - 5S = 0.462 
mPa ■ s = 462 fj .Pa • s. From Eq. (9-12.2) 


1 - 433.2/536.78 
1 - 370.93/536.78 


0.624 


With a = 0.1175, Eq. (9-12.1) gives 

In t, = ln(0.1175 X 462) 


ln(462) 


ln(0.1175 X 462) 

v = 185 p.Pa • s = 0.185 cP 

„ 0.185 - 0.188 

Error = - X 100 = -1.6% 

0.188 


A more general estimation method would logically involve the extension of the 
high-pressure gas viscosity correlations described in Sec. 9-6 into the liquid region. 
Two techniques have, in fact, been rather widely tested and found reasonably ac¬ 
curate for reduced temperatures above about 0.5. These methods are those of 
Chung, et al. (1988) and Brule and Starling (1984). Both methods use Eq. (9-6.16), 
but they have slightly different coefficients to compute some of the parameters. The 
Chung, et al. form is preferable for simple molecules and will treat polar as well 
as nonpolar compounds. The Brule and Starling relation was developed primarily 
for complex hydrocarbons, and the authors report their predictions are within 10% 
of experimental values in the majority of cases. The Chung, et al. method has a 
similar accuracy for most nonpolar compounds, but significantly higher errors can 
occur with polar, halogenated, or high-molecular weight compounds. In both cases, 
one needs accurate liquid density data, and the reliability of the methods decreases 
significantly for T r less than about 0.5. The liquids need not be saturated; subcooled 
compressed liquid states simply reflect a higher liquid density. The Chung, et al. 
technique was illustrated for dense gas ammonia in Example 9-12. The procedure 
is identical when applied to high-temperature liquids. 


Discussion 

The quantity of accurate liquid viscosity data at temperatures much above the nor¬ 
mal boiling point is not large. In addition, to test estimation methods such as those 
of Chung, et al. or Brule and Starling, one needs accurate liquid density data under 
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the same conditions which apply to the viscosity data. This matching makes it 
somewhat difficult to test the methods with many compounds. However, Brule and 
Starling developed their technique so that they would be coupled to a separate 
computation program using a modified BWR equation of state to provide densities. 
They report relatively low errors, and this fact appears to confirm the general ap¬ 
proach (See also Brule and Starling, 1984). Hwang, et al. (1982) have proposed 
viscosity (as well as density and surface tension) correlations for coal liquids. 

Regardless of what high-temperature estimation method is chosen, there is the 
problem of joining both high- and low-temperature estimated viscosities should that 
be necessary. 


9-13 LIQUID MIXTURE VISCOSITY 


Essentially all correlations for liquid mixture viscosity refer to solutions of liquids 
below or only slightly above their normal boiling points; i.e., they are restricted to 
reduced temperatures (of the pure components) below about 0.7. The bulk of the 
discussion below is limited to that temperature range. At the end of the section, 
however, we suggest approximate methods to treat high-pressure, high-temperature 
liquid mixture viscosity. 

At temperatures below T r ~ 0.7, liquid viscosities are very sensitive to the 
structure of the constituent molecules (See Sec. 9-11). This generality is also true 
for liquid mixtures, and even mild association effects between components can often 
significantly affect the viscosity. For a mixture of liquids, the shape of the curve 
of viscosity as a function of composition can be nearly linear for so-called ideal 
mixtures. But systems that contain alcohols and/or water often exhibit a maximum 
or a minimum and sometimes both (Irving, 1977a). 

Almost all methods to estimate or correlate liquid mixture viscosities assume 
that values of the viscosities of the pure components are available. Thus the methods 
are interpolative. Nevertheless, there is no agreement on the best way to carry out 
the interpolation. Irving (1977) surveyed more than 50 equations for binary liquid 
viscosities and classified them by type. He points out that only very few do not 
have some adjustable constant that must be determined from experimental mixture 
data and the few that do not require such a parameter are applicable only to systems 
of similar components with comparable viscosities. In a companion report from the 
National Engineering Laboratory, Irving (1977a) has also evaluated 25 of the more 
promising equations with experimental data from the literature. He recommends 
the one-constant Grunberg-Nissan (1949) equation [see Eq. (9-13.1)] as being 
widely applicable yet reasonably accurate except for aqueous solutions. This NEL 
report is also an excellent source of viscosity data tabulated from the literature. 
Other data and literature sources for data may be found in Aasen et al. (1990), 
Aucejo, et al. (1995), supplementary material of Cao, et al. (1993), Franjo, et al. 
(1995), Kouris and Panayiotou (1989), Krishnan, et al. (1995, 1995a), Kumagai and 
Takahashi (1995), Petrino, et al. (1995), Stephan and Hildwein (1987), Stephan and 
Heckenberger (1988), Teja, et al., (1985), and Wu, et al. (1998). 


Method of Grunberg and Nissan (1949) 

In this procedure, the low-temperature liquid viscosity for mixtures is given as 

j n n 

In Tj m = S In 17 , + ~ 2 S XiXjG v (9-13.1) 

i ^ i=l j—i 

or, for a binary of 1 and 2, 
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In rj m = x 1 In t) 1 + x 2 In tj 2 + x 1 x 2 G l2 (9-13.2) 

since G„ = 0. In Eqs. (9-13.1) and (9-13.2), x is the liquid mole fraction and G„ 
is an interaction parameter which is a function of the components i and j as well 
as the temperature (and, in some cases, the composition). This relation has probably 
been more extensively examined than any other liquid mixture viscosity correlation. 
Isdale (1979) presents the results of a very detailed testing using more than 2000 
experimental mixture datum points. When the interaction parameter was regressed 
from experimental data, nonassociated mixtures and many mixtures containing al¬ 
cohols, carboxylic acids, and ketones were fitted satisfactorily. The overall root 
mean square deviation for the mixtures tested was 1.6%. More recently, Isdale, et 
al. (1985) proposed a group contribution method to estimate the binary interaction 
parameter G, ; at 298 K. 

The procedure to be followed is: 

1. For a binary of i and j, select i by following the priority rules below. (j then 
becomes the second component.) 

a. i = an alcohol, if present 

b. i = an acid, if present 

c. i = the component with the most carbon atoms 

d. i = the component with the most hydrogen atoms 

e. i = the component with the most —CH 3 groups 
Gjj = 0 if none of these rules establish a priority. 

2. Once the decision has been made which component is i and which is j, 
calculate 2A,- and 2A■ from the group contributions in Table 9-12. 

3. Determine the parameter W. (If either i or j contains atoms other than carbon 
and hydrogen, set W = 0 and go to step 4.) Let the number of carbon atoms in i 
be Nj and that in j be N r 

(0.3161)(AT. - N) 2 

w = -- N ' + N J ~ (0.1188)(Af - Al.) (9-13.3) 

4. Calculate G, ; from 

G y = 2A,- - 2A,- + W (9-13.4) 

Gy is sometimes a function of temperature. However, existing data suggest that, 
for alkane-alkane solutions or for mixtures of an associated component with an 
unassociated one, Gy is independent of temperature. However, for mixtures of non¬ 
associated compounds (but not of only alkanes) or for mixtures of associating com¬ 
pounds, Gy is a mild function of temperature. Isdale, et al. (1985) suggest for these 
latter two cases, 


Gy(T) = 1 - [1 


Gy (298)] 


573 - T 
275 


(9-13.5) 


where T is in kelvins. 

Example 9-21 Estimate the viscosity of a mixture of acetic acid and acetone at 323 
K (50°C) that contains 70 mole percent acetic acid. Isdale, et al. quote the experimental 
value to be 0.587 cP, and, at 50°C, the viscosities of pure acetic acid and acetone are 
0.798 and 0.241 cP, respectively. 


Downloaded from Digital Engineering Library @ McGraw-Hill (www.digitalengineeringlibrary.com) 
Copyright © 2004 The McGraw-Hill Companies. All rights reserved. 

Any use is subject to the Terms of Use as given at the website. 





VISCOSITY 


VISCOSITY 


9.79 


TABLE 9-12 Group Contributions for G, at 298 K 


Group 

Notes 

Value of A,. 

—ch 3 


-0.100 

>CH, 


0.096 

>CH— 


0.204 

>C< 


0.433 

Benzene ring 
Substitutions: 


0.766 

Ortho 


0.174 

Meta 


— 

Para 


0.154 

Cyclohexane ring 


0.887 

—OH 

Methanol 

0.887 


Ethanol 

-0.023 


Higher aliphatic alcohols 

-0.443 

V 

n 

II 

o 

Ketones 

1.046 

—Cl 


0.653-0.161 N a 

—Br 


-0.116 

—COOH 

Acid with: 



Nonassociated liquids 

-0.411 + 0.06074A C 


Ketones 

1.130 


Formic acid with ketones 

0.167 


/V (:1 = number of chlorine atoms in the molecule. 

N c = total number of carbon atoms in both compounds. 


solution First we must estimate G, at 298 K. Component i is acetic acid (priority rule 
b). Since the mixture contains atoms other than carbon and hydrogen (i.e., oxygen), W 
= 0. Then, with Table 9-12, 

2A,. (acetic acid) = —CH 3 + —COOH = -0.100 + 1.130 = 1.030 
2A, (acetone) = (2)(—CH 3 ) + >C=0 = (2)(-0.100)+ 1.046 = 0.846 

With Eq. (9-13.4), 

G, = 1.030 - 0.846 = 0.184 at 298 K 
At 50°C = 323 K, we need to adjust G, with Eq. (9-13.5). 

0,(323 K) = 1 - a - 0.184X573 - 323) = „ 25g 

Then, using Eq. (9-13.2), 

In 7) m = (0.7) In (0.798) + (0.3) In (0.241) + (0.7)(0.3)(0.258) = -0.531 
r) m = 0.588 cP 

This estimated value is essentially identical with the experimental result of 0.587 cP. 

To summarize the Isdale modification of the Grunberg-Nissan equation, for each 
possible binary pair in the mixture, first decide which component is to be labeled 
i and which j by the use of the priority rules. Determine EA,. and EA- by using 
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Table 9-12 and W from Eq. (9-13.3), if necessary. Use Eq. (9-13.4) to calculate G i} . 
Correct for temperatures other than 298 K, if necessary, with Eq. (9-13.5). With 
the values of Gy so determined, use either Eq. (9-13.1) or (9-13.2) to determine the 
viscosity of the liquid mixture. This technique yields quite acceptable estimates of 
low-temperature liquid mixture viscosities for many systems, but Table 9-12 does 
not allow one to treat many types of compounds. Also, the method does not cover 
aqueous mixtures. 


UNIFAC-VISCO Method (Chevalier, et al., 1988; Gaston-Bonhomme, et ah, 
1994) 

Gaston-Bonhomme, Petrino and Chevalier have modified the UNIFAC activity co¬ 
efficient method (described in Chap. 8) to predict viscosities. In this method, vis¬ 
cosity is calculated by 


A*? £c 

In r) m = 2 .r, ln(r,,y ; ) - In V m + -±r + (9-13.6) 

The combinatorial term is the same as in the UNIQUAC model (see Table 8-8) and 
is calculated by 


A *g EC v , <h , z v , 0, 

p T = z x t In + > q,x, In 

i\.i i Xj 2 i 

(9-13.7) 

where z is the coordination number, equal to 10, 0, and </>, are the molecular surface 
area fraction and molecular volume fraction, respectively, given by 

x ; q ; 

0, = 

2j Xjqj 
j 

(9-13.8) 

and 


XT- 

4, = v 

2, -Vi 

j 

(9-13.9) 

where q t , the van der Waals' surface area, and r,, the van der Waals’ volume of 
component i, are found by summation of the corresponding group contributions. 
Thus, if rtf’ is the number of groups of type k in the molecule i. 

q t = X n \kQk 

k 

(9-13.10) 

r i = S n?R k 

(9-13.11) 


k 


where Q k and R k are the constants representing the group surface and size and are 
given in Table 9-13. These values match the UNIFAC values in Table 8-23 in cases 
where groups are defined the same. The residual term in Eq. (9-13.6) is calculated 
by 
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TABLE 9-13 UNIFAC-VISCO, Group 
Volume and Surface Area Parameters 


Group k 

R k 

Q k 

CH 2 , CH 2cy 

0.6744 

0.540 

ch 3 

0.9011 

0.848 

CH^ 

0.5313 

0.400 

Cl 

0.7910 

0.724 

CO 

0.7713 

0.640 

coo 

1.0020 

0.880 

OH 

1.0000 

1.200 

CHjOH 

1.4311 

1.432 


where 


and 


&*g ER 

RT 


2 x i In 


7i 


(9-13.12) 


In T,* ff = 2 n< kl In yf ~ In y * (0 ] (9-13.13) 

k 


In yf 


(\ In 1 V ^ 'XT/' * 1 

1 

g 

©' 

U 111 m * mk J 

^ 2 e„ 

n — 


2 X k Q k 


(9-13.14) 

(9-13.15) 


In Eq. (9-13.15), @ m is the surface area fraction in the mixture of groups and X m 
is the mole fraction in the mixture of groups. Except for the minus sign in Eq. 
(9-13.12), these last four equations are identical to those in the UNIFAC method 
described in Chap 8. However, the groups are chosen differently and the interaction 
parameters are different and are calculated by 


= exp ( - (9-13.16) 

Values of a nm are given in Table 9-14. y'f is the activity coefficient of group k in 
a mixture of groups in the actual mixture, and ■yjf® is the activity coefficient of 
group k in a mixture of groups formed from the groups in pure component i. Groups 
in branched hydrocarbons and substituted cyclic and aromatic hydrocarbons are 
chosen as follows 
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Type of compound 

Actual group 

Representation 

branched cyclic 

>CH—CH 3 

2 CH, groups 


>CH cy —CH 3 

1 CH, cy + 1 CH 2 


>C cy —(CH 3 ) 2 

1 CH acy + 2 CH, 

aromatic 

>c„ —ch 3 

1 CH,„ + 1 CH, 


Table 9-15 compares results calculated with the UNIFAC-VISCO method to ex¬ 
perimental values. Of all the methods evaluated, the UNIFAC-VISCO method was 
the only one that demonstrated any success in predicting viscosities of mixtures of 
compounds with large size differences. The method has also been successfully 
applied to ternary and quaternary alkane systems. The average absolute deviation 
for 13 ternary alkane systems was 2.6%, while for four quaternary systems it was 
3.6%. The method is illustrated in Example 9-22. 


Example 9-22 Use the UNIFAC-VISCO method to estimate the viscosity of a mixture 
of 35.4 mole% «-decane (1) and 64.6 mole% n-tetratetracontane, C^H,,, (2) at 397.49 
K. The experimental viscosity and density (Aasen, et al., 1990) are 3.278 cP and 0.7447 
g/cm 3 . 


solution From Aasen, et al. (1990), t/i = 0.2938 cP, 17 , = 4.937 cP, V l = 220 
g/cm 3 , and V 2 = 815.5 g/cm 3 . 


V„ 


2 x t M t 

Pm 


0.354 X 142.28 + 0.646 x 619.16 
0.7447 


604.7 cm 3 /mol 


In decane, there are 8 CH 2 groups and 2 CH 3 groups. In tetratetracontane, there are 42 
CH 2 groups and 2 CH 3 groups. Equations (9-13.10) and (9-13.11) give 


TABLE 9-15 Comparison of Calculated and Experimental Liquid Mixture Viscosities 


1 st 

component 

2 nd 

component 

Al 

T, K 

l?exp 

mPa-s 

Ref* 

'Hcalc 

mPa-s 

% 

deviation 

"-CioH 22 

^■^60^122 

0.749 

384.1 

3.075 

1 

2.309 

-25 



0.749 

446.4 

1.423 

1 

1.275 

-10 

72-Ci 0 H 2 2 


0.354 

368.8 

5.286 

1 

5.256 

- 0.6 



0.354 

464.1 

1.465 

1 

1.654 

13 



0.695 

374.1 

2.318 

1 

1.960 

-15 

butane 

squalane 

0.839 

293.1 

1.060 

2 

0.8812 

-17 

ethanol 

benzene 

0.5113 

298.1 

0.681 

3 

0.6403 

- 6.0 

acetone 

benzene 

0.3321 

298.1 

0.4599 

4 

0.4553 

- 1.0 

acetone 

ethanol 

0.3472 

298.1 

0.5133 

5 

0.4860 

-5.3 


*References: 1, Aasen, et al. (1990); 2. Kumagai and Takahashi (1995); 3. Kouris and Panayiotou (1989); 
4, Petrino, et al. (1995); 5, Wei, et al. (1985) 
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r, = 8 x 0.6744 + 2 X 0.9011 = 7.1974 
r 2 = 42 x 0.6744 + 2 X 0.9011 = 30.127 


q t = 8 X 0.54 + 2 X 0.848 = 6.016 
q 2 = 42 X 0.54 + 2 x 0.848 = 24.376 
Equations (9-13.8) and (9-13.9) give 




0.354 X 6.016 

0.354 X 6.016 + 0.646 X 24.376 


0.1191, 


ft, = 0.8809 


<t>i 


_0.354 X 7.1974_ 

0.354 X 7.1974 + 0.646 x 30.127 


0.1158, 


= 0.8842 


Equation (9-13.7) is used to calculate the combinatorial contribution 


A*g £c 

RT 


0.354 In 


0.1158 

0.354 


0.646 In 


0.8842 

0.646 


+ 5 


0.354 X 6.016 In 


0.1191 

0.1158 


+ 0.646 X 24.376 In 

0.8842 


= 0.1880 


In the mixture of groups, with CH, designated by subscript 1 and CH 3 by subscript 2: 


X. 


8 x 0.354 + 42 X 0.646 

8 X 0.354 + 42 x 0.646 + 2 X 0.354 + 2 X 0.646 


0.9374, 


X 2 = 0.0626 


Equation (9-13.15) gives 


0 , 


_0.9374 X 0,54_ 

0.9374 x 0.54 + 0.0626 x 0.848 


0.9051, 


© 2 = 0.0949 


Equation (9-13.16) gives 



66.53 \ 
298 ) 


0.7999, 


'Pfj = exp 


/ 709.5 \ 
\ 298 ) 


10.81 


Equation (9-13.14) gives 


In yf = 0.54 


1 - ln(0.9051 + 0.0949 x 10.81) — 


0.9051 


0.9051 

0.0949 X 0.7999 
0.9051 x 0.7999 + 0.0949 


0.0949 x 10.81 


= -0.1185 


Similarly. In y* = -3.3791 
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In the mixture of groups from pure component 1, X\'\ = 0.8, and = 0.2. Using 
Eq. (9-13.15) for pure component 1 then gives 


00 ) = 


In-yfM = 0.54 


0.8 X 0.54 


0.8 x 0.54 + 0.2 X 0.848 
1 - ln(0.7181 + 0.2819 X 10.81) - 


= 0.7181, O' 1 ' = 0.2819 

0.7181 


0.7181 + 0.2819 X 10.81 

0.2819 X 0.7999 
0.7181 x 0.7999 + 0.2819 


= -0.4212 


Similarly, •y) !<1) = -1.0479. In pure component 2, the results are @f’ = 0.9304, 0t 2) 
= 0.0696, In yf <2 > = -0.07655, and In y< 2 > = -4.1201. Equation (9-13.13) gives 

In y* K = 8 x (-0.1185 + 0.4212) + 2(-3.3791 + 1.0479) = -2.241 

In y* R = 42(—0.1185 + 0.07655) + 2(—3.3791 + 4.1201) = -0.2799 

Finally, the residual contribution is calculated with Eq. (9-13.12) 

A * o ER 

—= -(-0.354 x 2.241 - 0.646 x 0.2799) = 0.9741 


Equation (9-13.6) is now used to calculate the mixture viscosity 

In t)„, = 0.354 In (0.2938) + 0.646 In (4.937) + 0.354 In (220) 
+ 0.646 In (815.5) - In (604.7) - 0.1880 + 0.9741 

U„ = 3.385 cP 


Error = 


3.385 - 3.278 
3.278 


X 100 = 3.26% 


Method of Teja and Rice (1981, 1981a) 

Based on a corresponding-states treatment for mixture compressibility factors (Teja, 
1980; Teja and Sandler, 1980) (See chap. 5), Teja and Rice proposed an analogous 
form for liquid mixture viscosity. 

co — 

In(i 7 m e,„) = ln(i 7 £) (R1) + Lln(Tj£) (R2 ' - ln(i 7 £) (Rn ] ^ (R2) _ ^ (R1) (9-13.17) 

where the superscripts (Rl) and (R2) refer to two reference fluids, rj is the viscosity, 
co the acentric factor, and e is a parameter similar to £ in Eq. (9-4.15) but defined 
here as: 


y2/3 

7 C _ 

C t l mY' 2 


(9-13.18) 


The variable of composition is introduced in four places: the definitions of co,,,, V cm , 
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T cm , and M m . The rules suggested by the authors to compute these mixture para¬ 
meters are: 


V = 


2 2 x iXjV ci 

i j 

2 2 x i x jT C ijV C ij 


rri ■' 

^ cm t t 

v cm 

M m = 2 xMi 

i 

“■ = 2 -w 

i 

T j (VI? + O 3 

8 

7 ..y.. = ill ( /■ T V V )' 

cij cij tij \ ci cj ’ ci cj / 


(9-13.19) 

(9-13.20) 

(9-13.21) 

(9-13.22) 

(9-13.23) 

(9-13.24) 


i [ijj is an interaction parameter of order unity which must be found from experi¬ 
mental data. 

It is important to note that, in the use of Eq. (9-13.17) for a given mixture at a 
specified temperature, the viscosity values for the two reference fluids 17"' 11 and 
17 (R2) are to be obtained not at T, but at a temperature equal to /~|(7' r ) ,RI) /7’ m ] for 
(Rl) and T[(T c ) (R2) /T cm ] for (R2). T cm is given by Eq. (9-13.20). 

Whereas the reference fluids (Rl) and (R2) may be chosen as different from the 
actual components in the mixture, it is normally advantageous to select them from 
the principal components in the mixture. In fact, for a binary of 1 and 2, if (Rl) 
is selected as component 1 and (R2) as component 2, then, by virtue of Eq. (9- 
13.22), Eq. (9-13.17) simplifies to 


ln(i 7 ,„E,„) = Xi ln(i7e ) 1 + * 2 ln(i7s ) 2 


(9-13.25) 


but, as noted above, is to be evaluated at T(T cl /T cm ) and tj 2 at T(T c2 /T cm ). 

Our further discussion of this method will be essentially limited to Eq. (9-13.25), 
since that is the form most often used for binary liquid mixtures and, by this choice, 
one is assured that the relation gives correct results when x 1 = 0 or 1.0. In addition, 
the assumption is made that the interaction parameter 1 //- is not a function of tem¬ 
perature or composition. 

The authors claim good results for many mixtures ranging from strictly nonpolar 
to highly polar aqueous-organic systems. For nonpolar mixtures, errors averaged 
about 1%. For nonpolar-polar and polar-polar mixtures, the average rose to about 
2.5%, whereas for systems containing water, an average error of about 9% was 
reported. 

In comparison with the Grunberg-Nissan correlation [Eq. (9-13.1)], with G- 
found by regressing data, Teja and Rice show that about the same accuracy is 
achieved for both methods for nonpolar-nonpolar and nonpolar-polar systems, but 
their technique was significantly more accurate for polar-polar mixtures, and par¬ 
ticularly for aqueous solutions for which Gmnberg and Nissan’s form should not 
be used. 
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Example 9-23 Estimate the viscosity of a liquid mixture of water and 1,4-dioxane at 
60°C when the mole fraction water is 0.83. For this very nonideal solution, Teja and 
Rice suggest an interaction parameter i/r- = 1.37. This value was determined by re¬ 
gressing data at 20°C. 

solution From Appendix A, for water, T c = 674.14 K, V c = 55.95 cm 3 /mol, and 
M = 18.02; for 1,4-dioxane, T c = 587 K, V c = 238 cm 3 /mol, and M = 88.11. Let 1 
be water and 2 be 1,4-dioxane. With Eq. (9-13.18), e, = (55.95) 2 ' 3 /[(647.14)(18.02)] 1 ' 2 
= 0.135; e 2 = 0.169. From Eq. (9-13.19), 

„ „ [(55.95) 1 ' 3 + (238) 1 ' 3 ] 3 

V cm = (0.830) 2 (55.95) + (0.170) 2 (238) + (2)(0.830)(0.170) x -^-—2-—- 

8 

= 80.93 cm 3 /mol 
and with Eq. (9-13.20), 

T cm = {(0.830) 2 (647.14)(55.95) + (0.170) 2 (587)(238) + (2)(0.830)(0.170)(1.37) 
[(647.14)(55.95)(587)(238)] 1/2 } /81.29 = 697.9 K 
M m = (0.830)(18.02) + (0.170X88.11) = 29.94 
So, with Eq. (9-13.18), 


£ m 


(80.93) 273 

[(697.9)(29.94)] 1/2 


0.129 


Next, we need to know the viscosity of water not at 333.2 K (60°C), but at a temperature 
of (333.2X647.14)/697.9 = 309.0 K (35.8°C). This value is 0.712 cP (Irving, 1977a). 
[Note that, at 60°C, 17 (water) = 0.468 cP.] For 1,4-dioxane, the reference temperature 
is (333.2)(587)/697.9 = 280.3 K (7.1°C), and at that temperature, r, = 1.63 cP (Irving, 
1977a). Again this value is quite different from the viscosity of 1,4-dioxane at 60°C, 
which is 0.715 cP. Finally, with Eq. (9-13.25), 

ln[(r),„)(0.129)] = (0.830) ln[(0.712)(0.135)] + (0.170) ln[( 1.63)(0.169)] 


= -2.163 
r) m = 0.891 cP 

The experimental viscosity is 0.89 cP. 

Although the agreement between the experimental and estimated viscosity in Ex¬ 
ample 9-23 is excellent, in other composition ranges, higher errors occur. In Fig. 9-14, 
we have plotted the estimated and experimental values of the mixture viscosity over 
the entire range of composition. From a mole fraction water of about 0.8 (weight 
fraction = 0.45) to unity, the method provides an excellent fit to experimental results. 
At smaller concentrations of water, the technique overpredicts rj m . Still, for such a 
nonideal aqueous mixture, the general fit should be considered good. 


Discussion 

Three methods have been introduced to estimate the viscosity of liquid mixtures: 
the Grunberg-Nissan relation [Eq. (9-13.1)], the UNIFAC-VISCO method [Eq. 
(9-13.6] and the Teja-Rice form [Eq. (9-13.24)]. The Grunberg-Nissan and Teja- 
Rice forms contain one adjustable parameter per binary pair in the mixture. The 
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0 0.2 0.4 0.6 0.8 1.0 

Mole fraction water 

FIGURE 9-14 Viscosity of water and 1,4-dioxane at 333 K. 
Line is Eq. (9-13.25) with >p tJ = 1.37; o experimental. ( Irving, 
1977a) 


UNIFAC-VISCO method is predictive, but limited in the types of compounds to 
which it can be applied. The method correctly predicts the behavior of the 
methanol-toluene system which demonstrates both a maximum and minimum in 
the viscosity vs. concentration curve (Hammond, et al., 1958). An approximate 
technique is available to estimate the Grunberg-Nissan parameter G, ; as a function 
of temperature [Eq. (9-13.5)] for many types of systems. Teja and Rice suggest that 
their parameter i/^ is independent of temperature-at least over reasonable temper¬ 
ature ranges. This latter technique seems better for highly polar systems, especially 
if water is one of the components, and it has also been applied to undefined mixtures 
of coal liquids (Teja, et al., 1985; Thurner, 1984) with the introduction of reference 
components [See Eq. (9-13.16)]. The UNIFAC-VISCO method has been success¬ 
fully applied to ternary and quaternary alkane mixtures (Chevalier, et al., 1988) but 
otherwise, evaluation of the above methods for multicomponent mixtures has been 
limited. 

The above three methods are by no means a complete list of available methods. 
For example, Twu (1985. 1986) presents an equation to estimate the viscosity of 
petroleum fractions based on the specific gravity and boiling point. This method is 
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particularly useful for cases where the exact chemical composition of a mixture is 
unknown. Allan and Teja (1991) have also presented a method applicable to petro¬ 
leum fractions and (Chhabra, 1992) presents a method for mixtures of hydrocar¬ 
bons. Chhabr and Sridhar (1989) extend Eq. (9-11.6) to mixtures. For the treatment 
of electrolyte solutions, the reader is referred to Lencka, et al. (1998). Cao, et al. 
(1993) presented a UNIFAC-based method but our testing did not reproduce their 
excellent results in a number of cases. Other mixture correlations are reviewed in 
Monnery, et al. (1995) as well as the 4 lh edition of this book. For an example of 
gases dissolved in liquids under pressure, see Tilly, et al. (1994). 

An equation developed by McAllister (1960) has been used successfully to cor¬ 
relate data for binary as well as multicomponent mixtures (Aminabhavi, et al., 1982; 
Aucejo, et al, 1995; Dizechi and Marschall, 1982a; Noda, et al.. 1982). For binaries, 
the McAllister (1960) equation has been written to contain either two or three 
adjustable parameters. For ternary mixtures, the equation has been used with one 
(Dizechi and Marschall, 1982a) or three (Noda, et al., 1982) ternary parameters in 
addition to the binary parameters. Dizechi and Marschall (1982) have extended the 
equation to mixtures containing alcohols and water and Asfour, et al. (1991) have 
developed a method to estimate the parameters in the McAllister equation from 
pure component properties. Because of the variable number of parameters that can 
be introduced into the McAllister equation, it has had considerable success in the 
correlation of mixture viscosity behavior. 

Lee, et al. (1999) used an equation of state method to successfully correlate the 
behavior of both binary and multicomponent mixtures. Nonaqueous mixtures re¬ 
quired one parameter per binary while aqueous mixtures required two parameters 
per binary. One drawback of their method is the non-symmetrical mixing rule used 
for multicomponent aqueous mixtures (Michelsen and Kistenmacher, 1990). The 
equation of state structure allowed the method to be successfully applied to liquid 
mixtures at high pressure. 

To finish this section, we again reiterate that the methods proposed should be 
limited to situations in which the reduced temperatures of the components com¬ 
prising the mixture are less than about 0.7, although the exact temperature range 
of the Teja-Rice procedure is as yet undefined. 

Should one desire the viscosity of liquid mixtures at high pressures and tem¬ 
peratures, it is possible to employ the Chung, et al. (1988) method described in 
Sec. 9-7 to estimate high-pressure gas mixture viscosities. This recommendation is 
tempered by the fact that such a procedure has been only slightly tested, and usually 
with rather simple systems where experimental data exist. 


Recommendations to Estimate the Viscosities of Liquid Mixtures 

To estimate low-temperature liquid mixture viscosities, either the Grunberg-Nissan 
equation [Eq. (9-13.1) or (9-13.2)], the UNIFAC-VISCO method LEq. (9-13.6)] or 
the Teja-Rice relation [Eq. (9-13.17) or (9-13.25)] may be used. The Grunberg- 
Nissan and Teja-Rice methods require some experimental data to establish the value 
of an interaction parameter specific for each binary pair in the mixture. In the 
absence of experimental data, the UNIFAC-VISCO method is recommended if 
group interaction parameters are available. The UNIFAC-VISCO method is partic¬ 
ularly recommended for mixtures in which the components vary greatly in size. It 
is possible to estimate the Grunberg-Nissan interaction parameter G tj by a group 
contribution technique and this technique can be applied to more compounds than 
can the UNIFAC-VISCO method. All three methods are essentially interpolative in 
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nature, so viscosities of the pure components comprising the mixture must be 
known (or in the Teja-Rice procedure, one may instead use reference fluids of 
similar structure rather than the actual mixture components). The errors to be ex¬ 
pected range from a few percent for nonpolar or slightly polar mixtures to 5 to 
10% for polar mixtures. With aqueous solutions, neither the Grunberg-Nissan form 
nor the UNIFAC-VISCO method is recommended. 


NOTATION 


a* group contribution sum; Eq. (9-4.22) 

b a excluded volume, (2/3 )t tN 0 ct 3 , Eq. (9-6.1) 

C„ heat capacity at constant volume, JV(mol-K); C,, structural contribu¬ 

tion in Eq. (9-4.22) and Table 9-3 C br , number of carbon atoms in 
a branch 

D diffusion coefficient, cm 2 /s or m 2 /s 

F c shape and polarity factor in Eq. (9-4.11); F° p , low-pressure polar cor¬ 

rection factor in Eq. (9-4.18); F° Q , low-pressure quantum correction 
factor in Eq. (9-4.19); F P , high-pressure polar correction factor in 
Eq. (9-6.10); F g , high-pressure quantum correction factor in Eq. 
(9-6.11) 

g :j radial distribution function, Eq. (9-7.14) 

Gj, G 2 parameters in Eqs. (9-6.21) and (9-6.22); G, ; parameter in Eq. (9-13.1) 
A *g EC combinatorial contribution to viscosity in Eq. (9-13.6) 

A *g ER residual contribution to viscosity in Eq. (9-13.6) 

k Boltzmann’s constant 

L mean free path 

m mass of molecule 

M molecular weight 

n number density of molecules; number of components in a mixture 

N number of carbon atoms or parameter in Eq. (9-11.2); A N, structural 

contribution in Eq. (9-11.5) and Table 9-10; N 0 , Avogadro’s number 
P pressure, N/m 2 or bar (unless otherwise specified); P c , critical pres¬ 

sure; P r , reduced pressure, PIP ( : P vp , vapor pressure; A P r = 

(P - P vp )/P c 

q ; surface area parameter for molecule i 

Q polar parameter in Eq. (9-6.5); Q k , surface area parameter of group k 

r distance of separation; r ; , volume of molecule i 

R gas constant, usually 8.314 J/(mol ■ K), 7? , parameter in Eq. (9-6.4); 

R k , volume parameter of group k 
T* kT/s 

T temperature, K; T c , critical temperature; T r , reduced temperature, 

7/7); T b , boiling point temperature; T fp , melting point temperature 
v molecular velocity 

V volume, cm 3 /mol; V c , critical volume; V r , reduced volume, V/V c , or 

in Eq. (9-6.3), V/V 0 \ V,„ hard packed volume used in Eq. (9-6.3) or 
parameter in Eq. (9-11.6) 
x mole fraction, liquid 

y mole fraction, vapor; parameter in Eq. (9-6.20) 
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Y parameter in Eq. (9-6.9) 

Z compressibility factor; Z c , critical compressibility factor; Z,. Z 2 , param¬ 

eters in Eqs. (9-6.6) to (9-6.8) 

Greek 

y orientation factor in the Brule-Starling method. Table 9-8, or obtain 

from Brule-Starling (1984) 

A correction term in Eq. (9-6.17) 

e energy-potential parameter; variable defined in Eq. (9-13.18) 

77 viscosity (usually in micropoises for gas and in centipoises for liq¬ 

uids); 17 °, denotes value at low-pressure (about 1 bar); r \ c , at the 
critical point; 77 °, at the critical temperature but at about 1 bar; 17 *, 
17 **, parameters in Eqs. (9-6.19) and (9-6.23), 77,, reduced viscos¬ 
ity, defined in either Eq. (9-4.13) or (9-6.4), r/ ; ,, at the normal boil¬ 
ing point temperature 

6„ surface area fraction of molecule i or group k 

k polar correction factor in Eq. (9-4.11), see Table 9-1 

A thermal conductivity, W/ (m • K) 

77, dipole moment, debyes; 7 7 ,., dimensionless dipole moment defined in 

either Eq. (9-4.12) or Eq. (9-4.17) 
v kinematic viscosity, 17 /p, nr/s 

£ inverse viscosity, defined in Eq. (9-4.14) or Eq. (9-4.15); f r , inverse 

viscosity defined in Eq. (9-6.13) 
fa volume fraction of molecule i 

p density (usually mol/cm 3 ); p c critical density; p r , reduced density, 

p!Pc 

cr molecular diameter, A 

(// intermolecular potential energy as a function of r 

’T radial distribution function 

faj interaction parameter in Eq. (9-13.24) 

co acentric factor, Sec. 2-3 

collision integral for viscosity 

Subscripts 

i, j, k components i, j, k 

1 , 2 components 1 , 2 

L liquid 

171 mixture 

SL saturated liquid 
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